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ABSTRACT 


This report presents the results of work on communications 
systems modeling and covers three different areas of modeling. 

The first of these deals with the modeling of signals in communi- 
cation systems in the frequency domain and the calculation of 
spectra for various modulations. These techniques are applied 
in determining the frequency spectra produced by a unified carrier 
system, the down- link portion of the Command and Communications 
System (CCS). 

The second modeling area covers the modeling of portions of 
a communication system on a block basis. A detailed analysis 
and modeling effort based on control theory is presented along 
with its application to modeling of the automatic frequency con- 
trol system of an FM transmitter. Both linear and nonlinear 
models for the system are developed. 

A third topic discussed is a method for approximate model- 
ing of stiff systems using state variable techniques. Such sys- 
tems are characterized by having one or more very fast time con- 
stants along with much slower time constants. A technique for 
removing "fast” roots associated with fast time constants from the 
state equations characterizing a linear system is described, and 
an example presented. 

A number of the computer routines developed and used in the 
above efforts are listed in the report Appendices. 
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I. INTRODUCTION 


This report presents the results of a portion of the technical 
program conducted under Contract NAS8-20054, and covers three dif- 
ferent areas of communication system modeling. These are: (1) the 

modeling of signals in communication systems in the frequency domain 
and the calculation of spectra for various modulations, (2) the 
modeling of portions of a communication system on a block basis 
using control theory techniques, and (3) a technique for approxi- 
mate modeling of a "stiff" linear system that can reduce the solu- 
tion time significantly and still produce acceptable estimates of 
the system response. 

Section II presents results of work in signal modeling and 
calculation of spectra, in particular, determination of frequency 
spectra produced by a unified carrier system, the down-link por- 
tion of the Copimand and Communications System (CCS). Both cal- 

-JT> 

culated spectra and.iexperimentally measured spectra are presented 
for comparison. Also discussed is a technique for modeling the 
transfer function of filters. 

Section III describes a detailed analysis and modeling effort 
based on control theory and its application to modeling of the 
automatic frequency control system of an FM transmitter. The 
modeling was done on a block basis and both linear and nonlinear 
models for the system are developed. The linear model was verified 
with an analytic solution; nonlinear effects were introduced one 
at a time by remodeling a single block. The models were used to 
investigate various AFC characteristics, such as stability, tran- 
sient time, and frequency sensitivity. Results of these investi- 
gations are presented. 

Section IV discusses the method for approximate modeling of 
stiff systems using state variable techniques. Such systems are 
characterized by having one or more very fast time constants along 
with much slower time constants. Straightforward digital simula- 
tion of such systems is generally unsatisfactory in that excessive 
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round-off errors occur. A technique for removing "fast” roots 
associated with fast time constants from the state equations is 
described, and an example presented. 

Since the three sections deal with three different topics 
which are not closely related, the conclusions applicable to each 
effort are included in each individual section. 

A number of the computer routines developed and used in the 
above efforts are listed in the Appendices. All of the programs 
are written in FORTRAN-V for use on a UNIVAC-1108 computer. 
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II. SOME SPECIAL CONSIDERATIONS IN MODELING 
UNIFIED CARRIER SYSTEMS 


A, Introduction 

An investigation was made of the spectra of both baseband and 
RF signals in the CCS transponder when the modulating signals are 
a PCM data stream, a biphase modulated telemetry subcarrier, or 
both. The specific objective of the investigation was to develop 
techniques applicable to unified carrier systems which would per- 
mit ready determination of the effects on radiated spectra of 
such factors as number and types of modulating signals, modulation 
indices, and baseband filtering of signals. The work included 
modeling of various signals, modeling of filter transfer functions, 
and automating procedures for calculating and plotting spectra. 
During the course of the work one particular application was made; 
evaluation of the effect on the operation of the CCS transponder 
of replacing the PN range code on the down-link with a wide band 
data signal. This task was* performed on a quick-response basis and 
has been previously reported [l]. This report covers the particular 
areas investigated and the techniques developed on a more general 
basis. Recommendations are included on additional tasks that need 
to be carried out. 

Particular aspects of the CCS have been investigated in two 
previous studies. The work reported here makes use of techniques 
developed on these programs. These include mathematical derivation 
of expressions for certain baseband and RF spectra, development of 
computer programs for computing the spectra, automated plot routines 
for constructing graphs of spectra and of time waveforms from com- 
puted results, and an efficient mechanization of the fast Fourier 
transform (FFT) . For detailed discussions of these developments, 
reference should be made to the reports of these studies [2,3]. 

Some of the more pertinent points applicable to this study are 
reviewed briefly in the following sections. 


3 



B, The CCS Down-Link Configuration of Interest 

The CCS is a unified carrier system which combines the func- 
tions of command, communication, and ranging on a single S-band 
carrier for each direction of transmission [4]. The up-link portion 
of the system usually carries both command and ranging signals and 
was the subject of an earlier study [2]. The modulation applied 
to the down- link is more complex. The transponder aboard the space- 
craft demodulates the up- link signal and re-modulates it onto the 
down-link. In addition, a biphase modulated 1024 kHz telemetry 
subcarrier is also modulated onto the down-link. The down-link 
portion of the system was also investigated in a previous study 
[3]. More detailed discussions of the system and its characteris- 
tics are given in the reports on these earlier efforts [2,3]. 

In the CCS transponder, the telemetry subcarrier is normally 
biphase modulated by data at a rate of 72 kilobits per second. In 
the study reported here a data rate of 27.5 kilobits per second 
was assumed to correspond to a proposed use of the transponder. 

Also this study does not consider an actual turned-around up- link 
signal, but substitutes a PCM signal having approximately the same 
clock rate as the range code in lieu thereof. Thus for purposes 
of this study, the signals applied to the S-band modulator of the 
down link consisted of (1) a PCM signal replacing the turned-around 
up-link signal, and (2) a 1024 kHz subcarrier which was biphase 
modulated with a data rate of 27.5 kilobits per second. Of special 
interest in this study were the time waveforms and spectra of the 
original baseband signals, the transfer function of various filters, 
time waveforms and spectra of the filtered baseband signals, and 
the S-band spectra produced by these signals. 

The techniques developed here should be useful in the general 
area of communication system modeling and evaluation. 
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C, Theoretical Considerations and Computer Programs 

This section presents a brief discussion of the theoretical 
basis for calculating frequency spectra and time waveforms. In 
addition, the calculation of selected filter transfer functions 
is discussed. Computer programs developed for the calculations 
are also discussed, and listings of these programs are included in 
the Appendix. For more detailed discussion of the spectral cal- 
culation techniques, see Technical Reports Nos. 3 and 7 [2,3]. 


1. Baseband Spectral Calculations 

The input signal considered was a time waveform repre- 
senting the PCM signal, the modulated telemetry subcarrier, or 
their sum. The input signal is evaluated at equally spaced time 
intervals which gives the sampled version of the signal as a 
real-time function. The baseband frequency spectrum of the time 
waveform for the continuous signal case is given by 


x(f) 



x(t)e 


-2njf t 


dt 


( 1 ) 


where x(f) is the frequency domain representation, and x(t) is 
the time domain representation. For a sampled signal the discrete 
Fourier transform (DFT) must be used. An efficient means of com- 
puting the discrete Fourier transform (DFT) is provided by the FFT. 
The mechanization of the FFT used in these investigations is de- 
scribed in Appendix C of Technical Report No. 7 [3]. The frequency 

domain representation, A , of the samples, X, , of the time waveform 

r k. 

is given by 
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W = j^e — is the phase function of the DFT. Index "r" 
is referred to as the frequency of the transform component • The 
transform provides the spectral components, A^, corresponding 
to the series of time samples, X^. The inverse transform is 
obtained in a similar manner with the elimination of the 1/N 
term in front of the summation and by use of a positive sign in 
the exponential. 

The continuous signal representation of a filter transfer 
function is given by 

H(jcu) = ACaOe^^ = a + jb (4) 


where A(u)) is the amplitude response of the filter and 0(ui). is 
the phase function. The latter is given by 



(5) 


where a and b are the real and imaginary parts of H(juo). 

For the continuous case the filtered time function, g(t), may be 
obtained from the inverse transform of the product of H(joo) and 
F(jcu) as 


g(t) 



• F(jo)) 



( 6 ) 


by 


For the discrete case, the filtered time function is given 
n-1 

C (k) = ^ [G (i) • H(i) J W lk (7) 

i=o 


where C(k) is the sampled output time function, G(i) is the trans- 
form of the sampled input time function, and H(i) is the transform 

ik 

of the filter impulse response. W is the Fourier transform phase 
function. 


6 



This procedure is easily mechanized for computation. The 
filtered time function can be generated by transforming the 
sampled time function to the frequency domain, evaluating the fil- 
ter transfer function at these transform frequencies, producing 
the product of the transform frequency component with the value 
of the transfer function (in general both complex numbers) , and 
re transforming the result to the time domain. 


2. Radio Frequency Spectral Calculations 

The time waveform representing a phase modulated carrier 
such as that present in the CCS can be expressed as 


e 

pm 


(t) 



e 


-Jtt) «: -jPt(t) 

+ e e 


m 


( 8 ) 


where 

A = the carrier peak amplitude, 
od^ = the carrier radian frequency, 

(3 = the phase modulation index, and 
\|r(t) = the modulation time function. 


The two terms in brackets in Equation (8) represent the positive 
and negative frequency components of the modulated carrier, respec- 
tively. The carrier frequency in the CCS is high enough so that 
only the positive frequency terms need to be considered. Also, 
since only relative magnitudes are of interest in the spectral cal- 
culations, the time function representing the modulated carrier 
can be considered to be 

j(JU t jp\|t(t) 

e!_ (t) = e e . (9) 

pm 

The Fourier transform of the modulated carrier signal is 


E'(^) = 6 (to — uu ) * F(cu) 
c 


( 10 ) 
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where 6 (id - a> c ) represents a delta function at the carrier frequency, 
the symbol denotes convolution, and F(oj) = The 

Fourier transform, F(o)), of the modulating signal can be obtained 
by evaluating the exponential function by knowing the modulation 
indices and the values of the time function. Evaluating this expo- 
nentialfunction and taking the FFT of the result yields the radio 
frequency spectrum. (Subroutine FOFT which evaluates is 

given in Appendix A. ) 

3. Filter Transfer Functions 

Several filter transfer functions were modeled for the 
wideband PCM data down-link studies. Two representations of filter 
transfer functions were used. One of these expresses the filter 
transfer characteristics as a function of the ratio of the frequency 
at which the transfer function is to be evaluated to the filter 
cutoff frequency (3 dB point). For this representation, component 
values of the filter are not required to evaluate the transfer 
function -- only the cutoff frequency of the filter need be con- 
sidered. This method of generating the transfer function was used 
in modeling a RC low-pass filter and is discussed in Technical 
Report No. 7 [3]. 

The other representation of a transfer function considers the 
actual component values of the filter. The loop equation of the 
filters are written and solved in determinate form. Examples of 
the application of this technique to several filters are also dis- 
cussed in Technical Report No. 7 [3]. The application of this 
technique to the 1024 kHz bandpass telemetry filter in the Gener- 
alized Concept Receiver will be discussed here. This filter has 
a transformer coupled stage and presents a good example for demon- 
strating computer modeling of a representative filter. 

The schematic diagram of the telemetry bandpass filter is 
shown in Figure 1. The equivalent circuit of the filter in which 
the input transformer is replaced by a ,, T" network is shown in 
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Figure 2a. The equivalent network of Figure 2b is used for writing 
the loop equations. Each impedance shown in Figure 2b was repre- 
sented by its appropriate combination of circuit elements as shown 
in Figure 2a. The various impedances are 
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Writing the loop equations and solving for the transfer func- 
tion gives 


H 



( 12 ) 
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The expressions (11) relating the impedances with numbered sub- 
scripts to circuit parameters include resistances (e.g., Rll, R13, 
etc.) which determine the Q of the inductances. The Q used for 
all inductances had a value of 100. Also included is the self 
capacitance of the larger inductances. These appear in (11) as 
the C values with the double subscripts. The self capacitance 
of the coils was experimentally obtained using the relationship 



where C = self capacitance of the coil, 

u 

- capacitance required to resonate the coil at a 
frequency f^, and 

C 2 = capacitance required to resonate the coil at a 
frequency f^ - 2 f . 

The self capacitance of the coils was determined experimentally 
from the filter inductances to be 

C n = 20.7 pf, 

C 22 = 4#27 pf > 

C 33 = 16.07 pf, 

= 4.1 pf , and 

C 55 = 16.93 pf. 

The self capacitance was included in the circuit only for coils 
with large inductances for which the self capacitance produced 
self resonance in the 5 MHz region. The self resonance of the 
smaller coils was in the region of 20 MHz, well removed from 
the filter passband. 

The telemetry filter transfer function was evaluated using 
subroutine TLMFLT included in Appendix B. The results of the 
computations are shown in Figure 3 which displays the normalized 
amplitude response as well as the phase response. The amplitude 
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AMPLITUDE (OB) , PHASE ( OEGREES ) 



Figure 3. Amplitude and phase response of COS telemetry- 
bandpass filter. 



response of the filter obtained in the laboratory is shown in 
Figure 4. This display of the amplitude response of the filter 
was obtained using a spectrum analyzer and a constant amplitude 
signal generator whose frequency was swept very slowly. The 
response at zero frequency in Figure 4 is the "zero beat" of the 
spectrum analyzer and does not represent the value of the trans- 
fer function of the filter which is zero at this frequency. The 
calculated response shows good agreement with the experimentally 
obtained value. 

The effect of the telemetry bandpass filter transfer func- 
tion on the square wave modulated telemetry subcarrier is shown 
in Figures 5 through 7. Figure 5 shows the computed biphase 
modulated telemetry subcarrier modulated at the nominal CCS data 
rate of 72 kilobits per second. Figure 6 shows the frequency 
spectrum of the signal. The telemetry filter transfer function 
was then applied to this spectrum and the results retransferred 
to the time domain. Figure 7 shows the time waveform at the out- 
put of the filter. 

D. Discussion of Binary PCM Baseband Systems 

Binary PCM data waveforms can take several forms depending 
on the method used to indicate a one or a zero. The TRIG tele- 
metry standards specify seven different methods of bit coding. 

Two common methods of bit coding were considered in this analysis. 
These are specified in the standards as NRZ-Level, in which a one 
is represented by one level and a zero by the other, and Biphase- 
Level, in which a one is represented by a one-zero transition 
during a bit time and a zero is represented by a zero-one transi- 
tion during a bit time. NRZ-Level will be referred to here simply 
as non-return- to-zero (NRZ) and Biphase-Level as return- to-zero (RZ) . 

For NRZ data having a bit time of one microsecond the bit 
rate is one megabit per second and the fundamental clock rate is 
500 kHz. For RZ data with the same time for each state of the 
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( a ) HORIZONTAL CALIBRATION 200 kHz PER DIVISION. 
VERTICAL CALIBRATION 10 dB PER DIVISION. 



( b ) HORIZONTAL CALIBRATION 0.5 MHz PER DIVISION 
STARTING AT 1 MHz. 

VERTICAL CALIBRATION 10 dB PER DIVISION. 

Figure b. Experimentally determined amplitude response of 
the CCS bandpass telemetry filter. 
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Figure 






signal one bit time would occupy two microseconds and the bit 
rate would be one half that of the NRZ mode. NRZ data with an 
alternating string of ones and zeros would have a fundamental 
frequency of 500 kHz. Similarly, a string of ones or zeros with 
RZ data would also produce a square wave at 500 kHz. If the data 
were a perfect square wave passed through a linear filter, har- 
monic components of the data would exist only at odd harmonics 
of the data fundamental frequency. For an alternating bit pat- 
tern for RZ data the fundamental frequency would be one-half of 
that for NRZ data and again harmonics would exist at odd har- 
monics of the fundamental. For strings of ones or zeros in RZ 
data it would appear as NRZ data with an alternating bit pattern 
since a transition is forced to occur in each time slot assigned 
to a bit. Data such as those just discussed would have no detri- 
mental effect on the telemetry system operating in a bandwidth 
of approximately 150 kHz centered at 1024 kHz. Although NRZ 
data would provide a higher data rate, RZ data is often preferred 
because of the additional signal transitions forced into the data 
leading to better performance of tape recorders and bit synchronizers. 

E. Examples of Some Calculations and Experimental Measurements 

for a Proposed use of a Unified Carrier System 

The techniques developed for the time waveform and frequency 
spectrum modeling of communication systems, such as the CCS down- 
link, were used for the analysis of the signals involved in a pro- 
jected use of a unified carrier system such as the CCS. Parts of 
the results of this investigation were reported earlier in a special 
technical memorandum which provided, on a quick- look basis, data 
relative to the problem [l]. Some additions to the data avail- 
able at that time have been made. This intended use of a unified 
carrier system provides a good example of the modeling of systems 
using frequency domain techniques. 

The investigation of the spectral content of baseband and 
radio frequency signals in the CCS transponder down-link for spe- 
cific baseband signals was of interest. 


20 


These signals were a biphase modulated 1024 kHz telemetry sub- 
carrier signal and a 500 kilobit per second split phase modulated 
PCM data signal. For the case of interest the modulation on the 
telemetry subcarrier was nominally a 13.75 kHz square wave simulat- 
ing NRZ data at a rate of 27.5 kilobits per second. The split 
phase modulated PCM signal was simulated by a 500 kHz square wave. 
This simulation represents a continuous string of ones or zeros 
in the split phase modulated data stream. A bit time in this sig- 
nal occupies two microseconds. For a split phase signal, a one 
in the input data stream is represented by a one-to-minus-one 
transition at the midpoint of the bit interval, while a zero is 
represented by a minus-one-to-one transition at the midpoint of 
the bit interval. Thus, a continuous stream of ones or zeros 
results in a square wave having a period of two microseconds or 
a frequency of 500 kHz. 

Calculation of the spectra of the various signals involved 
was accomplished using programs developed earlier for CCS spectral 
studies. To adapt the input signals to the FFT, slight adjust- 
ments were made in the frequencies of some of the signal s» The 
biphase modulating signal frequency was adjusted to be 13.65333 kHz 
and the wideband PCM signal frequency adjusted to 498.3467 kHz. 

This made all signals periodic on a period of two cycles of the 
biphase subcarrier modulating signal (146.484375 p^sec) with 73 
cycles of the PCM signal and 150 cycles of the 1024 kHz subcarrier 
frequency present during this period. These adjustments were made 
to avoid some unknowns about the discrete transform when the trans- 
form coefficients do not match the frequencies present in the in- 
put signal. 

Measured spectra were obtained from laboratory measurements 
for comparison with computed spectra. The modulation indices 
were set experimentally by determining the peak-to-peak input 
voltage to the phase modulator of a sine wave signal at either 
500 kHz or 1024 kHz required to give the desired index. The 
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phase modulation index for the two modulating signals was then 
set by matching the peak-to-peak amplitude of each of the signals 
to the amplitude of the appropriate sine wave. 

The experimental test setup is shown in Figure 8. The tele- 
metry subcarrier portion of the baseband signal was generated by 
applying the biphase modulating signal and the telemetry subcarrier 
signal to the proper inputs of a balanced modulator. The output 
signal of the balanced modulator is the desired biphase modulated 
telemetry subcarrier signal. The 500 kilobit per second data 
signal was obtained from a pulse generator adjusted to give a 
500 kHz square wave signal. These two signals were added in an 
oscilloscope and the vertical output of the oscilloscope was used 
as the input to the CCS transponder phase modulator. The ampli- 
tudes of the two signals could be adjusted independently by use 
of the controls on the two channel vertical amplifier of the 
oscilloscope . 

Modulation indices used in the calculated and experimental 

data were 0.7 radian for the wideband PCM signal (Pp C ^) and 1*0 

radian for the telemetry subcarrier signal (|3 ). Figure 9 shows 

s c 

the time waveform of such a modulating signal, and Figure 10 
shows its frequency spectrum. Figure 11 shows the same time 
waveform after filtering through two ideal RC low pass filters 
having a cutoff frequency of 1 MHz each. The spectrum of the 
output of the two low-pass sections is shown in Figure 12. This 
signal serves as the input to the down- link phase modulator. 

Experimentally obtained baseband spectra are shown in Fig- 
ure 13. Figure 13a shows the biphase modulated subcarrier spec- 
trum for the above example; Figure 13b is that of the wideband 
PCM signal; and Figure 13c shows the sum of the two spectra. 

An example of the time waveform for a random selection of 
the 73 bits possible in the PCM signal is shown in Figure 14. 

The selection of the bit pattern was made using a pseudo-random 
number generator in a Univac 1108 computer. That is, each of 
the 73 bits of the PCM signal present in the time waveform was 
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(a) BIPHASE MODULATED (SQUARE 
WAVE MODULATION) 1024 kHz 
TELEMETRY SUBCARRIER. 

VERTICAL lOdB/DIV. 
HORIZONTAL 200 kHz/DIV. 



(b) 500 KILOBIT/SEC SPLIT PHASE 
SIGNAL CONSISTING OF ALL 
ONES OR ZEROS. 

(SIMULATED BY A 500 kHz 
SQUARE WAVE.) 

NOTE: Sidebands around square 
wave fundamental and 3rd 
harmonic frequencies were 
found to be produced by 
the square wave signal 
generator. 




Figure 13. Baseband spectra of modulating signals. Amplitude 

relationships for modulation indices of 3 =1.0 

and 6 = O.T. SC 

pcm 
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selected in a random fashion with a probability of 0*5. Figure 15 
shows the baseband spectra of the signal shown in Figure 14. Note 
the disappearance of the single line in the spectrum at 500 kHz, 
and the spreading of the spectrum around this frequency. Note 
also the null in the spectrum of the PCM data signal near the tele- 
metry subcarrier frequency. 

Next the S-band spectra of some of these signals were investi- 
gated. Figure 16 shows the calculated S-band spectrum for the bi- 
phase modulated telemetry subcarrier signal only. The calculated 
S-band spectrum for both modulating signals present is shown in 
Figure 17. Experimentally obtained S-band spectra for each modula- 
tion present separately and both at the same time are shown in 
Figure 18. 

Ideally, the effect of phase demodulation should be evaluated 
to complete the study, but time did not permit developing modeling 
techniques for the phase demodulator. For the particular applica- 
tion investigated, experimental assessment of the demodulator 
action was made, and the results are shown in Figure 19. This 
figure shows the output of the telemetry bandpass filter in the 
subcarrier phase demodulator in the Generalized Concept Receiver 
[5] for the case of each signal presented individually and both at 
the same time. Note that when both signals are present the second 
harmonic component of the 500 kilobit per second signal (imperfect 
square wave) is some 10 dB below the amplitude of the first side- 
band pair of the biphase modulated telemetry subcarrier signal. 

F, Conclusions 

Modeling studies have been carried out which were aimed at 
developing techniques and procedures applicable to unified carrier 
systems which would permit analytical determination of the effect 
of changing modulation types, modulating indices, etc. Techniques 
have been developed for modeling biphase modulated subcarrier sig- 
nals and PCM data streams. Automated procedures were developed 
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Figure l6. S-band spectrum of a biphase modulated telemetry 
subcarrier only, 6 = 1.0. No low-pass simulation 
of modulator characteristic. 



Figure 17. S-band spectrum of a biphase modulated telemetry subcarrier, 
6 = 1.0, and Manchester PCM signal consisting of all ones or 
zeros, B = 0.7 > no low-pass simulation of phase modulator 



(a) MODULATED (SQUARE WAVE) 
1024 kHz TELEMETRY SUB- 
CARRIER' 

VERTICAL lOdB/DIV. 
HORIZONTAL 360kHz/DIV. 



(b) 500 KILOBIT SPLIT PHASE 
PCM SIGNAL. 




3U 





(a) MODULATED 1024kHz 
TELEMETRY SUBCARRIER 
ONLY. 

VERTICAL lOdB/OIV. 
HORIZONTAL 200 kHz/OIV. 





Figure 19 • Baseband signal spectra at output of 102 L kHz 

Generalized Concept Receiver telemetry bandpass 

filter, 3 =1.0, and 6 =0.7. 

sc pcm 


35 



for plotting the time waveform of baseband signals; for computing 
both baseband and RF spectra; and for plotting the spectra. These 
procedures permit easy determination of the spectral content of 
a phase modulated carrier when both biphase modulated subcarriers 
and PCM data streams are used as modulating signals. 

A need exists for extending this work to modeling of a phase 
demodulator, so that the time waveform and spectral content of the 
recovered baseband signal can be predicted. Time available for 
this program did not permit modeling of the demodulator. 

During the course of the work a specific application of the 
techniques was made in examining the effects of replacing the PN 
range code on the CCS down- link with a wideband PCM signal. It 
was found that the effects of the PCM signal on the recovered 
telemetry waveform were small. 
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III. AUTOMATIC FREQUENCY CONTROL SYSTEM 
ANALYSIS AND SIMULATION 


A. Introduction 

An Automatic Frequency Control (AFC) loop for a typical 
direct Frequency Modulation (FM) transmitter is considered in 
this section. A digital simulation (nonlinear) of the AFC sys- 
tem is developed and, together with a linear analysis, is used 
to determine the relationships between system behavior and the 
parameters and other characteristics of the subsystems. Using 
such a simulation, the design engineer can evaluate and improve 
his preliminary system design without setting up the associated 
hardware and test equipment. Two methods of frequency sampling 
are considered. 

After a brief description of a direct FM transmitter and 
its AFC system, math models are described for the subsystems of 
the AFC loop. Using linearized models, important AFC system 
features such as stability, transient time, and frequency sensi- 
tivity are estimated analytically. These linear system results 
provide insight into the significance of subsystem parameters 
and serve as a test case for verifying the programming of the 
simulation . 

After the linear version of the simulation has been verified, 
certain linear blocks or subsystems in the simulation are replaced 
by nonlinear models. The Voltage Controlled Oscillator (VCO) in- 
put-output characteristic is modeled by a nonlinear curve-fit, 
and a limiter is included in the feedback path to simulate satura- 
tion in the error-processing electronics. 

The results of the analysis and simulation show how such 
parameters as the system gains, sampling period and duration, 
and the RC time constant of the sampling network influence the 
AFC system behavior. A numerical example is included. 
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B. System Description 


1, Direct FM Transmitters 

A typical direct FM transmitter with its associated 
Automatic Frequency Control (AFC) loop is shown in block diagram 
form in Figure 20 and again in more detail in Figure 21. In 
direct frequency modulation, the modulating voltage (input signal) 
varies the natural frequency of an LC oscillator whose nominal 
unmodulated frequency is the carrier frequency or some submultiple 
of it. The variation in natural frequency of the oscillator 
(ideally, proportional to the input signal) results when a change 
in the input voltage varies the capacitance of a voltage-dependent 
capacitor in the oscillator resonant circuit. 

This technique can produce large-index frequency modulation 
directly (i.e., without requiring large frequency multiplication 
ratios to produce the desired modulation index at the carrier fre- 
quency as in the Armstrong indirect frequency modulation technique), 
but since it uses an LC Voltage Controlled Oscillator (VCO) which 
tends to have drift in its nominal frequency, some method of 
stabilizing the carrier frequency to meet its stability require- 
ments is required. (A typical stability requirement may be 
0.005% at S-band.) To obtain such stabilities, the frequency 
of the LC oscillator is ‘'slaved 1 * to the frequency of a crystal 
reference oscillator through comparison and control circuitry 
(AFC loop). To accomplish this, a sample of the output signal 
can be heterodyned against a reference signal generated by multi- 
plication of the frequency of a crystal oscillator. The difference- 
frequency between the transmitter carrier signal and the multi- 
plied reference signal is applied to a discriminator (usually 
operating in the 1 MHz region) to produce an error signal for 
correction of the VOC nominal frequency. A low-pass filter 
(or sample and hold circuit) prevents the tracking-out of useful 
modulation components while providing for the correction of the 
long term frequency of the VCO. Frequency stabilization techniques 


38 




39 


Figure 20. Block diagram of FM transmitter. 
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for direct frequency modulation may operate to hold the mean fre- 
quency of the transmitted signal within a specified tolerance or 
they may sample the modulating signal at some repetitive position 
of the time waveform which is then held to some specific frequency 
(Figure 21). An example of the latter system would be a sampled 
automatic frequency control system used with a television modulat- 
ing signal. 

2 . AFC Loop 

In the analysis and simulation of the AFC system, the 
variable of interest is the frequency of the output signal rather 
than the output voltage waveform. Thus any RF block in Figures 
20 or 21 which does not influence the fundamental frequency of 
the waveform can be ignored in this study, regardless of any volt- 
age gain or waveform distortion the block (subsystem) might intro- 
duce. In the following, only that portion of the direct FM trans- 
mitter in the AFC loop will be discussed, and the discussion will 
treat frequency, rather than voltage, as the variable of interest. 

The heart of a direct FM transmitter is the voltage controlled 
oscillator (VCO) . Ideally, the VCO output frequency would be pro- 
portional to the modulating voltage; in reality it is a nonlinear 
function of the applied voltage. The VCO output (modulated signal) 
is amplified and perhaps frequency multiplied by succeeding stages 
in the transmitter. Although some of these stages included in 
the AFC loop may introduce phase shift and distortion in the mod- 
ulated signal, it is assumed that no change in the fundamental 
frequency is introduced (other than intentional frequency multi- 
plication) . Thus for simulation and analysis of the AFC system 
these stages may be modeled as constant gains. 

It is assumed that periodically in the input modulating sig- 
nal there is a reference level in the time waveform such as one 
level of a frequency shift keyed signal or the signal level during 
the sync pulse interval of a standard television video signal. 
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During the time of occurrence of this reference level , the output 
frequency of the transmitter is measured and compared with the 
desired nominal frequency. The difference frequency is then used 
as an error signal to adjust the bias voltage on the VCO so as to 
reduce the error in its output frequency. By sampling the error 
signal pulse during the center of its duration, the effects of 
the dynamics of the error detection system upon the AFC loop are 
minimized. Since the error signal is available only as a sampled 
signal and since a continuous bias for the VCO is required, there 
must be a sample-and-hold device between the error frequency sig- 
nal and the VCO bias input. Note that if the feedback (including 
sample and hold) contains no integration, the steady state error 
of the AFC system will be non-zero for a constant disturbance and 
will be inversely proportional to the AFC loop gain. The inclusion 
of an integrator would convert the AFC system to a control system 
which would reduce the steady-state error in the output frequency 
to zero for a constant disturbance caused by a frequency shift 
internal to the VCO. 

In the next section of this report, math models will be estab- 
lished for each block of the AFC loop. These models are used for 
both analysis and digital simulation of an automatic frequency con- 
trol system of a typical direct FM transmitter. 

C. Math Models 

In the following, math models for the individual blocks of 
the AFC loop shown in Figure 21 are described. Linearized models 
will be used in both the analysis and digital simulation of the AFC 
system, while the more complete nonlinear models will be used in 
the simulation only. 

1. Subsystem Models 

The Voltage Controlled Oscillator (VCO) produces a sinu- 
soidal voltage signal whose frequency is dependent upon the input 
voltage, V, in a nonlinear manner as shown in Figure 22. 
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Figure 22. Typical VCO characteristics. 


Since the VCO frequency is the variable of interest, it is not 
necessary to model the VCO with circuit details. The change in 
VCO frequency for a change in input voltage is assumed to occur 
"rapidly" with respect to other time-constants in the AFC loop 
so that the dynamics of the VCO may be ignored. 

The nonlinear model sketched in Figure 22 can be used in the 
simulation either by (1) finding an analytic expression which 
approximates f^(V) sufficiently close, or (2) by interpolating 
between data points from a measured f - V curve. 

Alternately, as indicated in Figure 22, the VCO characteris- 
tics can be approximated for "small 1 ’ deviations in V from its 
nominal value by the following linear expression: 

f = f + k (V - V ) . 
o c V c 

The associated block diagram for a linearized VCO model is as given 
in Figure 23. 

Although the Amplifier (see Figure 21) following the VCO 
definitely influences the waveshape of the signal (amplitude and 
phase), it is assumed not to influence the frequency of the signal. 
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f 0 = f c + k v (V-V c ) 



Figure 23- Linearized VCO model. 

Since the AFC loop analysis and simulation is concerned with only 
the frequency of the transmitted signal, the amplifier will be 
modeled here by a unity transfer function® Also, the output sig- 
nal Frequency Multiplication and Power Amplification blocks can 
be ignored for the purpose here since they are outside the AFC 
loop® 

The low-power, f ixed-f requency Reference Oscillator , whose 
frequency is a sub-harmonic of the desired nominal transmitter 
frequency, is modeled as a cons tant-f requency source. The ref- 
erence signal is fed into a Frequency Multiplier circuit with a 
multiplication factor n where the desired nominal transmitter 
frequency is (n + 1) f . The frequency multiplier is modeled 
as a gain factor, n. The VCO signal, after amplification, and 
the frequency multiplied reference signal are fed into a Mixer 
whose output signal has frequency equal to the difference fre- 
quency, (f Q - nf R ) ■> The deviation of this frequency from f 
represents the "error" in the nominal VCO output signal frequency. 
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The Limiter-Discriminator in Figure 21 is time-shared in the 
AFC loop to develop a voltage which represents alternately the fre- 
quencies (f - nf ) and f . Since the two signals may not have 
o R R 

the same amplitude, any Amplitude dependence 11 of the limiter- 
discriminator produces outputs as if the two input signals were 
processed by separate discriminators with different gains. This 
amplitude dependence will be modeled by using two discriminators 
whose gains are not necessarily the same. Nonlinear models for 
the discriminators could be used, but linear models with non- 
identical gains should be sufficient here. The (linearized) dis- 
criminator output voltage, V^, can be expressed in terms of the 
frequency of the input signal (f), the nominal discriminator fre- 
quency (f Q ) , and the nominal slope of the limiter-discriminator 
characteristic curve as follows: 

V 1 - k 01 (£ - £ o ) • <14) 

The linearized block diagram for the discriminator (modeled as a 
pair) is shown in Figure 24. 



Figure 2h. Discriminator Model 
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The difference between the two discriminator output voltages 

(or between the single discriminator's output voltages at the two 

different times) is proportional to the "frequency error" of the 

VC0„ After amplification, the difference voltage is combined with 

the "set" voltage, V^ ag , to form the feedback error signal. The 

voltage is set initially such that the output frequency, f Q , 

is at or near the desired nominal frequency for the reference value 

of the input voltage, V (see Figure 21). Provision is made in 

m 

the model to limit the signal at this point for modeling saturation 
throughout the feedback loop. 

Since the error signal is assumed available only on a sampled 
basis, and since the input voltage to the VCO must be continuously 
present, some Sample and Hold device must be used to buffer the 
feedback error signal. Depending upon how the time-sharing dis-~ 
criminator is implemented, the frequency error sample may be either 
a flat-topped pulse whose amplitude represents the error at the 
time the pulse began, or it may be a time-varying pulse whose ampli- 
tude follows the frequency error during the duration of the pulse. 
The former case, which allows a "deadbeat response," will be assumed 
here, and the other case discussed later. 

For the constant amplitude pulse, the error feedback can be 
modeled as two electronic switches, a zero-order sample-and-hold, 
and a simple RC charging circuit (Figure 25) . 

The initial value of the sampled waveform is held as a con- 
stant input to the RC circuit. The switch is assumed to close 

for t seconds once each T seconds. 

o 

This sample-and-hold circuit is piecewise linear and has a 
continuous output voltage, e^. However, the output is used in 
such a way in the AFC system that it is its new value, after 
has reopened, that is of interest. Accordingly, the sample-and- 
hold circuit can be modeled by the linear, first order, difference 
equation derived in the following. 

Consider one cycle or period for the sample-and-hold circuit. 
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SAMPLE-AND-HOLD 



Figure 25- Error feedback sampling model. 


Let the switch close at t = 0 and reopen at t = t Q and then stay 
open until t = T* Transient analysis shows that, while the switch 
is closed, the output voltage is given by 


e Q (t) = e Q (0) + [e in (0) - e Q (0)][l - e t//RC ], 0 < t < t Q , 

(15) 

and after the switch opens, the output voltage remains fixed at 

-t /RC 

e Q (t) = e Q (0) + [e. n (0) - e Q (0)][l - e ° ], t Q < t < T. 

(16) 


Thus , 


e o m 

where 


*.»> + C 1 - e o (0) J - 

-t /RC 

i ° 

1 - e 


This is generalized to the n 


th 


sampling period as 


(17) 


e o [(N + 1) T] = e o (NT) + ^ [e_ n (NT) - e Q (NT)] , 
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which can be expressed as a difference equation without the factor 
T as follows: 


e (N + 1) = e (N) + C,[e. (N) - e (N)] . (18) 

o o 1 in o 

Thus the updated sampler output will be equal to the current value 
plus a term proportional to the difference between the present in- 
put and output voltages. Note that this acts like a zero-order, 
or box-car, sample-and-hold device only if RC « t^ « T; that 
is, if I. However, even for t Q < RC, the AFC system closed 

loop response can be made similar to that of a box-car hold cir- 
cuit by using a large AFC feedback loop gain. This is shown in 
the analysis section of the report. 

2 o AFC Loop Model 

Using the math models just described for the blocks of 
Figure 21, the block diagram is redrawn as shown in Figure 26. 

A digital simulation program can be written directly from 
the block diagram, allowing a programming block for each diagram 
block. The program can of course use either the linearized models 
as shown in Figure 26 or the more complete nonlinear models de- 
scribed earlier. However, in order to verify a simulation, a 
reasonable test case is required; since the linear model response 
can be predicted analytically, it will be used as the test case. 
Individual blocks of the simulation will be replaced with non- 
linear models after the simulation program has been verified. 

The analysis of the linearized system (Figure 26) is given in 
the following section. 

Do Linear System Analysis 

The linearized model (Figure 26) of the AFC loop is analyzed 
in this section to (1) provide insight into the relationships be- 
tween system parameters and the effectiveness of the frequency 
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DISCRIMINATOR 


THE DIFFERENCE EQUATION IS V 4 (N+1) = V 4 (N) + C 1 [Vj (N) - V 4 (N)]. . 
NOTE: V B , AS =V C -(V M ) REF 


Figure 26. AFC loop model. 








control system and to (2) provide data for verifying or de-bugging 
the digital simulation program to be developed. The simulation 
will be developed initially using the linearized model, will be 
de-bugged; and will then have some of its (linear) blocks replaced 
by more complete models for further system evaluation. 

Since this system is a sampled-data system, the Laplace trans- 
form technique is not applicable for systems analysis. A similar 
technique using the Z-transform is normally used for sampled-data 
system analysis, but its use requires that the order of the denom- 
inator of the closed loop transfer function exceed that of the 
numerator by at least two [6]. Thus the Z-transform is not appli- 
cable here. 

Fortunately, the system is sufficiently simple that one can 
learn a lot about the system from its closed loop difference equa- 
tion. From this equation, such items of interest as stability of 
the closed loop system, sensitivity of the output frequency to 
variations in system parameters, and the time constant of the 
transient response can be determined. 

1. Closed°Loop Difference Equation 

Referring to the linearized model of the AFC system 
(Figure 26) the only subsystem which contains dynamics is the 
sampler block; all other subsystems are modeled simply as com- 
binations of gains and summations. Thus, with the exception of 

'k 

the difference equation modeling the sampler, the system math 
model is algebraic. Using the sampler difference equation and 
the algebraic relationships given by the linearized AFC block 
diagram, the closed-loop difference equation can be derived as 
follows . 


* 

A difference equation is the discrete analogy of a differential 
equation • 
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Let V^(N) refer to the value of the "sample and hold" out- 
put after the sampler switch has closed and reopened for the 

time. Similarly, let V (N) be the reference value of the modulat- 

m 

ing voltage just before the switch is closed and reopened for the 
(N + l)*** 1 time. 

The algebraic equation for the VCO is as follows: 
f (N + 1) = f + k [V(N + 1) - V ] . (19) 

O C V c 

But , 

V (N +1) = V (N + 1) - V. (N + 1) , 

m 4 

so that 

f (N + 1) = f - k V + k [V (N + 1) - V. (N + 1)] . (20) 

o cvcvm 4 

The difference equation for the sample and hold device is 

V 4 (N + 1) = [1 - C x ] V 4 (N) + C : V 3 (N) . (21) 

Substituting, 

f (N + 1) = f - k V + k V (N + 1) - k [(1 - C,)V. (N) + C.V 0 (N)] . 

o cvcvm v 14 1 j 

(22) 

Now, writing V, (N) in terms of V (N) and f (N) , writing V„ (N) in 
4 mo j 

terms of f Q (N), and substituting into Equation (22) will complete 
the development of the closed loop difference equation. 

Rewrite Equation (19) for N rather than (N + 1): 

f (N) = f - k V + k V (N) - k V. (N) , (23) 

o cvcvm v 4 

from which 

k V, (N) = -f (N) + f - k V + k V (N) . (24) 

v 4 o cvcvm v 7 
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From the block diagram, 


V N) - V bias + k A V 5 ’ 


‘ V bias + k A [V 1 ’ V • 


* V biaa + k A [k Dl <£ o (B) ‘ n£ R ' V ' k B2 ( V f D ) ] ’ 


or 


V B) - (k A K D1 )£ o OO + V blas - f R k A (nk D1+ k^) - k A f D (k D1 - k M ) 


(25) 


The loop is closed by substituting Equations (24) and (25) into 
Equation (22) to obtain 


f (N + 1) + [C.k k k - 1 + C ,] f (N) = k V (N + 1) 

o lvADl lo vm 

+ ([C 1 - 1] k v ) V ffl (N) + E , (26) 


where 


E - C.f • 
1 c 


C k V - 
1 v c 


C-k V + C k k (n k + k ) f 
1 v bias 1 v A Dl D2 R 


+ C l k v k A (k Dl- V f D 


(27) 


Equation (26) with the constant E as defined in Equation (27) is 
the first-order difference equation which describes the behavior 
of the Automatic Frequency Control System. Although an analytic 
expression (a finite sum series) exists for the solution to such 
a difference equation [7], the form of the expression is such that 
it is difficult to use in other than numerical examples. A better 
understanding of the significance of the parameters of the system 
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can be had by examining the difference equation directly for 
stability, sensitivity, and transient time. 

2. Stability Regions 

From the closed loop difference equation. Equation (26), 
much information about the response of the linearized model of 
the AFC system can be determined without actually solving the dif- 
ference equation. For example, one can determine the set of gains 
for which the closed loop system is stable, the "speed' 1 of response 
of the system, the steady state response, and the system sensitivi- 
ties (i.e., how the steady state output frequency varies with per- 
turbations in the system elements). 

To determine the transient response characteristics of a 
linear system, the driving functions may be set to zero and the 
resulting homogeneous differential equations or difference equa- 
tions analyzed. The homogeneous portion of Equation (26) is 

f c (N + 1) + [C 1 (k + 1) - 1] f o (N) = 0, (28) 


where 




v D1 A 


Regardless of the initial condition, f(N ® 0), the system 
described by Equation (28) can be seen to be stable for 

l [C 1 (k + 1) - 1]| < 1 (29) 

since 

f (N + 1) = - [C (k + 1) - 1] f (N) ; (30) 

o 1 o 

i.e., the magnitude of f (N) will be decreasing with each increment 
in N. Furthermore, the response is stable and non-oscillatory if 
-1 < [C^(k +1) - l] < 0, since f q (N + 1) would have the same sign 
as f Q (N). Similarly, 
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(31) 


0 < [C x (k + 1) - 1] < 1 

corresponds to a stable but oscillatory response with f decreasing 
in magnitude but reversing signs at each increment. Note the special 
form of response (called dead beat ) which results when 

[C 1 (k + 1) - 1] = 0. (32) 

Under this condition, the AFC loop reaches its "steady state" 
response within one sampling period of the occurrence of a step 
disturbance . 

The gain conditions for each type response are given in 
Table I. 


TABLE I 


TRANSIENT RESPONSES OF LINEAR AFC SYSTEM 


Loop Gain 

Factor 

Transient Response 


(k) 

[1 - C x (k + 1)] 


k> [f - 

c i 

1] 

Less than -1 

unstable , alternating 
sign 

k = [f - 

c i 

1] 

-i 

square wave 

Cr - 0 

■= k < ( f ■ 0 

-1 to 0 

stable , alternating 

v. i y 



sign 


0 

0 

deadbeat 

-1 < k < ( 


0 to +1 

stable, monotonic 





k = -1 


+1 

constant 

k < -1 


Greater than +1 

unstable, monotonic 

Note: 

-t /RC 

i ° 

= 1 - e 
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Consider a numerical example for which 


3 

R = 20 x 10 ohms, 

C = 0.1 x 10 ^ farads, and 

t = 1.5 x 10 seconds, 
o 

Then 

( i - l) - 1333 , and 

O' 2667 • 

Figures 27 shows the previously described stability regions. 

Typical AFC systems operate in the monotonically stable region 
with loop gain, k, well below the deadbeat value. The use of a dead- 
beat system, which gives the M best M frequency regulation as far as 
response time is concerned, requires a high loop gain which can 
make the system very susceptible to noise. Deadbeat control also 
requires more stages of error signal amplification; large pertur- 
bations can produce saturation in the loop. 

The stable-but-alternating portion of Figure 27 requires 
larger loop gain than deadbeat; it has no advantages over the 
stable and monotonic operation. 

As to the applicability of this linear model, it should be 
noted that the linearization of the VCO required the assumption 
that frequency perturbations were small. Thus the model is not 
applicable when the frequency deviation, Af, gets "large." The 
linear model does non the less provide a basis for choosing the 
loop gain. Figure 28 shows the transient response in frequency 
error for a step change in VCO nominal frequency for various values 
of AFC loop gain. 
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Figure 28. Transient responses of linear AFC system for 
various loop gains. 







3, Transient Response Parameters 


The linear model can also be used to predict an approxi- 
mate time constant, t, for the closed-loop AFC system response when 
k is in the stable, monotonic region. Since the approximate time 
constant is associated with the envelope of the steps in Af, it 
is significant only if several (> 10) sampling intervals are re- 
quired for 907o of the frequency error to be removed, i.e., only for 
relatively small loop gains. 

Recall that the homogeneous closed-loop difference equation 
is 

f (N + 1) - Y, £ (N) (33) 

o 1 o 


where 

Y x - 1 - C 1 (k + 1) 


-t /RC 

1 - (k + 1) (1 - e ° ) . 


For the stable monotonic response, 


0 < < 1 . 

By looking at the progression 


f o (T) - Vl f o (0) , 
f o (2I) " Vi f o m - Y 2 ! f 0 «» , 

f o (3I) ' £ 0 <2T> ' Yj f o (0) , 

etc o , 

one can see that the general homogeneous response is given by 

f 0 (NT) = Y N f Q (0), 
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or 


(t/T) 

f Q (t) = (y) f Q (0), where t = NT. 


But the approximate time constant, t, is defined such that* 


f 0 (T) = (e *) f Q (0) . 


Since 


(T/T) 


-u 


f c (T) = (Y) f Q (0) = (e *) f Q (0) , 


(t/t) -1 

Y - e , and 


T = 


-T 

In (Y) • 


T = 


-T 


In [1 - (k + 1) (1 - e 


-t /RC 


)] 


(34) 


Also, if t < < RC, 
’ o 9 


T to 


In [1 - (k + 1) t /RC] ; 


(35) 


if further [(k + 1) t /RC] » 1, then 

o 


T to 


(k + 1) (t /RC) 


(36) 


Thus by using either Equation (34) or Equations (35) or (36), 
one can choose the AFC system open-loop gain to obtain the desired 
transient response time constant. 

For higher gain (if ever desirable) and stable but alternating 
response, the overshoot is found from the homogeneous closed loop 
difference equation to be . 


Note that the approximate time constant may not be an integer 
multiple of the sample period, T. 
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4. Sensitivity Analysis 


In the analysis up to this point, the homogeneous dif- 
ference equation has been used to investigate the system transient 
response. Now that it has been shown that for a wide range of sys- 
tem parameters the closed-loop AFC system is stable, it is reason- 
able to ask how well the AFC system does its job. That is, how 
much will the VCO frequency deviate from nominal for a given para- 
meter change in the AFC loop. Although the term "frequency sta- 
bility" is commonly used to describe the sensitivity of the VCO 
frequency to system changes, one must be careful not to confuse 
this with stability in the transient sense; i.e., stable versus 
unstable . 

The system frequency stability or sensitivity can be deter- 
mined from the AFC closed-loop difference equation without actually 
solving the difference equation. Assume that step-changes in all 
system parameters occur at t = 0 (changes are assumed small enough 
that the linear model remains applicable). It is known that the 
system will converge eventually to a steady-state or constant fre- 
quency if the loop gain has been chosen to satisfy the inequality 


-1 < k < 



The steady-state frequency is determined from the AFC closed-loop 
difference equation by equating f^(N) and f^(N + 1). (Note that 
this can only be done after one is assured that the closed-loop 
system is not unstable.) The closed-loop difference equation is 
repeated here. 


f (N + 1) + [C,(k + 1) -1] f (N) = k V (N + 1) 
o I o v m 


+ (k [C, - 1] V (N) + E 
v 1 J m 


(37) 


where E is a function of the AFC system parameters, assumed con- 
stant at their new values: 
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E - C l' f c- “v^c + V bias ) + k v k A [(nk Dl + V f R 

+ (k M- V £ d ] > ' (38) 

For steady-state analysis, one has 

f (N) = f (N + 1) = f 
o o o 

and 

V (N) = V (N + 1) = V . 
mm m 

The difference equation then reduces to the following algebraic 
equation: 


C- (k + 1) f = C-k V 4* E , 
1 o 1 v m ’ 


or 


✓ k x 

f + 

o \k + 1/ m 


c x (k + 1) 


(39) 


This is the equation for the steady-state VCO frequency. Note 
that is a multiplying factor of E and thus does not influence 
the steady-state VCO frequency. That is, the system sampling 
parameter 


C 1 = 1 - e 


-t /RC 
o 


influences only the transient response of the AFC system. Sub- 
stituting for E from Equation (38) , the steady-state frequency 
is written as follows. 
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f 

o 


fjv_y v . v . v Vf — V +f p 

Mc + lA m c bias J V k + 1 J c R L 


k k A (nk Dl + k D2 


k + 1 


-] 


+ f. 


r Wsr k D2 

D L 


k + 1 


i] 


(AO) 


Nominally, set (experimentally) such that when the reference 

input voltage is present, the output frequency has its desired value 
of 


< f o> 

nom 

Also, the nominal values of other parameters are as follows: 

£ c - (n + 1) £ R , 

f D ’ V a " d 

k V2 ” k Dl = k o * 

From this expression for the steady state nominal output frequency, 
it can be seen that system parameter perturbations influence the 
frequency output as follows. Let 

Af - shift in VCO frequency for fixed input voltage, 

Af D = shift in reference oscillator frequency, 

Afp = shift in discriminator center frequency, 

AV = shift in reference input voltage, 
m 

AV^ = shift in bias voltage, and 

Af = resultant shift in nominal VCO frequency. 
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Then the resultant output frequency error is given by 


k k _ k k.Cnk^-H- k_ 0 ) 

kc - / v s AT7 / v v - , / 1 \ KC , v A D1 D2 KC 

af o ( k + 1* iV m ( k + 1^ AV b + ( k + 1* Af c + k + 1 A£ t 


. k v k A (k Dl- k D2 ) . . 

+ — im — Af D ■ 


(41) 


(Note: k = k k.k„„) 

v A D1 


which is approximated as follows for k » 1 and k^ p» k^ 


A£ o « (f ) AV m - <-*) AV b + -f + (n + 1) A£ r + (" I>1 k __ D2 > Af D 


A>1 ' k - 


Dl 


(42) 


The terms of Equation (42) are considered one-by-one in the 
following. The first shows the prime function of and a reason for 
the existence of the AFC loop; any shift in the reference voltage 
level of the input signal will result (steady state) in an error 
in the output VCO frequency that is (l/k) times that for an open- 
loop system. The second term shows that variations in the input 
bias voltage have a steady-state effect similar to shifts in the 
input reference level. The third term is the most important term. 
It shows that VCO frequency shifts due to any internal (to the VCO) 
changes will also be corrected by the AFC loop such that the resul- 
tant f error will be reduced by a factor of (l/k). The fourth 
term points out that percentage perturbations in the reference 
frequency appear directly (un-attenuated) in the VCO output fre- 
quency (note that f ^ (n + 1) f ) . The remaining term is small 
since it is of second order: both (k^- k^) and Af^ are normally 

small . 
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5, Summary of Analysis 

The preceding AFC system analysis serves two purposes. 

The insight gained in the functioning of the system enables one 
to better choose system parameters in design or re-design efforts 
and evaluate (without hardware) proposed changes in the system 
gains, etc. The linear analysis also can serve as a checkcase in 
the development of a computer simulation of the AFC system. Note 
that the simulation, after initial debugging, will use more com- 
plete, nonlinear models for the subsystems and will thus provide 
more realistic estimates of system behavior, especially for large 
frequency deviations which may produce saturation in the real 
system. The next section of this report describes the develop- 
ment of the computer simulation from the system block diagram. 

E. Simulation- Linear 

1. Development of Simulation Program 

The actual development of an AFC system simulation pro- 
gram requires only a small step beyond the development of the math 
model indicated by the block diagram in Figure 26 and does not 
depend directly upon the analysis in the previous section (except 
for de-bugging of the resulting program) . 

The program is written so as to repeatedly "go around" the 
AFC loop evaluating the (algebraic) transfer functions and updating 
the variables representing the states and/or outputs of each block* 
A "memory" of the previous (i.e., previous time around the loop) 
value of variables is required only for the variables associated 
with the dynamic transfer function(s). The predominance of alge- 
braic rather than dynamic blocks in the AFC system model makes 
the simulation program especially simple. 

The heart of the program is the repeated evaluation of only 
eight equations. Referring again to Figure 26, and beginning at 
the feedback point, the following equations are written: 
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V = V - v. , 
m 4 

f = f + k (V - V ) , 
o c v c 

f i = f o - nf R ’ 

V 1 = k Dl (f i " V ’ 

V 2 = k D2 (f R ‘ f D ) » 

v 5 = V 1 - v 2 , 

V. = V, . + k V , 

3 bias A 5 

v 4 (n+i)= v 4 (n) + [v 3 (n) - v 4 (n)] c x . 


(43) 


The basic linear AFC system simulation program consists of 
(1) an initialization section, (2) a DO LOOP to evaluate repeatedly 
the above equations, and (3) instructions for printing the simula- 
tion results, A listing of this basic program, including explana- 
tory comments, is given in Appendix E. Note that with the aid of 
the comments it is easy to see where to insert any desired non- 
linear model for a block or subsystem. 

A sample output from the basic program is also given in 
Appendix E. The loop gain for this run is 100. Note that the 
input step in occurred at the 10^ sample interval and that the 

system had essentially recovered (error frequency down by 907 o ) 
th 

by the 40 step. 


2. Verification of Linear Simulation 

For those loop gains for which the transient response 
is stable and monotonic, the system time constant is given by 
Equation (34) as derived earlier, or in graphical form by Figure 
29. Similarly, the resultant system steady state frequency error 
for a shift in nominal VCO frequency is given by Equation (41) 
and sketched in Figure 30. To verify the simulation (linear), 
the response time constants and steady state errors of simulation 
runs for several loop gains will be compared with the theoretical 
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APPROXIMATE TIME CONSTANT 
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Figure 29- AFC system approximate time constant. 


= VCO FREQUENCY DRIFT, WITHOUT AFC LOOP 
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Figure 30. Steady-state frequency error 




or expected values. The data continues to be that for an AFC 
loop with sampler parameters such that the characteristic of 
the sampler is described by = -.00075. For the comparisons, 

loop gains of 10, 25, 50, and 100 will be used. 

From the simulation runs from which the transient responses 
plotted in Figure 31 were obtained, the simulation "final 11 values 
of frequency error (approximate since settling was not complete) 
shown in Table II were read. These errors are to be compared with 
an initial error of 1.4 x 10 Hz; i.e., the ordinate of Figure 30 
is to be multiplied by this initial error before comparing. The 
simulation final errors are seen to be essentially in agreement 
with the predicted values; the small differences are a result of 
ending the simulation at a finite time. 

Having verified the steady state performance of the simula- 
tion program, one next investigates the transient behavior. In 
using the plot of the system time constant versus gain, note that 
this time constant is for the decay from the initial error fre- 
quency to the steady-state frequency error. That is, 

Af(T) = (Af ) + e" 1 [(Af). _ (Af) ] (44) 

s s c s s 

o 

TABLE II 


COMPARISON OF THEORETICAL AND SIMULATED FREQUENCY ERRORS 


k 

(Af ) 

ss 



Simulation 

Predicted 

10 

127596 

127272 

25 

54358 

53846 

50 

28014 

27451 

100 

14388 

13861 

500 

2796 

2794 
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FREQUENCY ERROR 
(NORMALIZED) 



NUMBER OF SAMPLING INTERVALS 

(a) 

Figure 31. Step response of the VCO frequency in a sampling AFC 
loop for various values of loop gain in the linear 
simulation . 







-.2 
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NUMBER OF SAMPLING INTERVALS 


(g) 


Figure 31 (Continued). 


72 










FREQUENCY ERROR (NO 



2 


3 


5 


NUMBER OF SAMPLING INTERVALS 





where t = T, the system time constant, 

(Af) 55 steady state frequency error, 
s s 

(Af) t = initial frequency error from which the system is 
o recovering. 

To verify the simulation transient response, the above expres- 
sion for the frequency error at one time-constant, Af(T), is eval- 
uated for each of several values of loop gain. (The predicted 

value of (Af) can be used.) The calculated values of Af(T) for 
ss 

k= 10, 25, 50, and 100 are tabulated in the second column of 
Table III. Referring then to the individual simulation runs, the 
time at which Af has the value Af(T) is noted and the elapsed time 
since the start of the decay, normalized by dividing by the sample 
period, is recorded in Column 3 of the table. Finally, the pre- 
dicted value of t/T is calculated from Equation (34) for each gain 
and is recorded in the fourth column of Table III. Note that the 
differences between the simulation and predicted values of t/T 
are less than one (step) and result from the use of a continuous 
expression, Equation (34), to approximate the discrete system. 

TABLE III 


COMPARISON OF THEORETICAL AND SIMULATED TIME CONSTANTS 


k 


t/T 




Simulation 

Predicted 

10 

595482 

121 

120.8 

25 

549068 

51 

50.8 

50 

532384 

26 

25.7 

100 

523793 

13 

12.7 
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As with the steady-state comparisons, the transient comparisons 
verify the correctness of the simulation. Now that the total AFC 
system is known to be correctly simulated, individual subsystems, 
or blocks, will be modeled in more detail in the simulation. Un- 
like the linear simulation, the later version will give results 
not available analytically. Specifically, the effects on the sys- 
tem performance of a limiter in the feedback path or a nonlinear 
VCO characteristic will be determined by simulation. 

F. Simulation - Nonlinear 

The linear simulation results show that as the loop gain is 
increased, both the transient response time and the steady-state 
frequency error go down, with the response time even reaching one 
sampling period for the gain referred to as "deadbeat gain." How- 
ever, for reasonable sized frequency errors, these larger gains 
would result in very high voltages in the AFC loop, and the linear 
model of the AFC loop is not applicable for high voltages. It 
still may be desirable to use high gains so that the frequency 
error will be small and the small perturbations will disappear 
rapidly, if the system nonlinearities do not introduce stability 
problems at the higher signal levels; to investigate this area, 
the simulation will be modified to include the two predominant 
nonlinearities . 

1. Limiter 

The saturation of the feedback signal processing stages 

will be modeled as a single limiter located between the output 

of the error signal amplifier and the input to the sample-and-hold 

circuit, at a point after the summing of V, . with the error sig- 

bias 

nal. The large-signal nonlinearity of the oscillator (VCO) input- 
output curve can be modeled in the simulation by replacing the 
linear VCO model with either a nonlinear analytic function fitted 
to actual VCO characteristics, or by using a table-look-up and 
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interpolation scheme on data points from an actual VCO charac- 
teristic curve. 

The saturation effect will be modeled by an ideal limiter. 

As long as the input signal has a magnitude less than the limit, 

L, the limiter has no effect on the signal; if, however, the in- 
put signal magnitude exceeds L, the output has the sign of the 
input signal but has magnitude L. 

Simulation runs with the limiter in the program with several 
values of loop gain, k, have been made and plots of the error fre- 
quency versus time are given in Figure 32. These responses should 
be compared with those for the linear model which have been given 
in Figure 31. The responses show that the effect of the limiter 
is to slow the recovery process to a rate determined by the limiter, 
until the error signal becomes small enough to be within the linear 
region of the limiter. Thus the shape of the response curve depends 
upon the amplitude of the disturbance. Since the (stability) 
effect of the limiter is similar to that of reduced gain, and 
since the AFC system does not become unstable for decreased posi- 
tive gains, the addition of the limiter would not be expected to 
introduce stability problems. This is born out by the simulation 
results. On the other hand, the limiting voltage must be set high 
enough so that expected variations in the reference value of the 

modulating signal, (V ) -, do not bias the system so much that 

m ret 

the sampler cannot supply the voltage required to offset V and 

m 

keep the VCO frequency near nominal. 

Figure 32 shows that the effect of including the limiter in 
the AFC loop simulation is to limit the slope of the response 
curve. Once the error signal is reduced below the saturation 
level, the response rate is determined by the time constant 
derived earlier; while still in saturation, the response rate is 
independent of the feedback gain. Note, however, that the value 
of Af at which the system saturates depends upon the loop gain. 


78 



.0075 

1.4 X 10 6 Hz/VOLT 
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Figure 32. AFC loop simulation response with limiter 
set at 12 volts. 




Thus the advantages of short transient time and low steady 
state frequency error associated with the higher loop gains in 
the linear model still apply. The only change in performance 
is the limiting of the response slope when the error frequency 
is large. Of course, the idea of a "deadbeat gain" is only 
applicable for the small frequency errors. 

2. Nonlinear VCO 

The Voltage Controlled Oscillator (VCO) characteris- 
tic curve (output frequency versus input voltage) was shown in 
Figure 22 as nonlinear but has been approximated for analysis 
and simulation thus far by a linear model. The VCO characteris- 
tic shown in Figure 33, which has an exaggerated nonlinearity, 
has been simulated in the AFC simulation to study the effects of 
the nonlinearity on the AFC system response. A listing of the 
AFC simulation program including the nonlinear VCO model and 
saturation is given in Appendix F. A typical simulation run is 
included. To simulate a particular VCO characteristic, one 
would simply replace the nonlinear equation with a curve-fit of 
the desired VCO characteristic. 

Figure 34 shows the frequency-error transient response 
curves for the nonlinear VCO model. Note that the system re- 
mains stable and that the form of the response is changed very 
little by the inclusion of the nonlinearity. As would be expected 
from the sketch of the VCO nonlinearity (Figure 33) , the VCO 
"gain" decreases as the input voltage increases above its nominal 
value and the gain increases as V drops below its nominal value. 

At the nominal point, the VCO "gain" is that used in the linear 
simulation. In the simulation runs, the input reference modulat- 
ing voltage is stepped by one volt from V = 7 to V = 8 volts. 

In the linear model this produces an initial frequency error of 
1.4 MHz; for the nonlinear VCO of Figure 33 the frequency shift 
is only 0.8 MHz. Note that if the voltage shift had been of the 
opposite polarity the nonlinear VCO frequency shift would have 
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Figure 3^. AFC simulation transient response with nonlinear VCO 




been -2.0 MHz while that for the linear one would have been only 
-1.4 MHz. 

From the transient responses in Figure 34 the effective gain 
reduction due to the VCO nonlinearity can be seen to decrease 
the magnitude of the initial slope. The effect is most easily 
seen for k = 25, but is present for all gains. In the linear 
model, the transient responses (when smoothed) were exponential 
decay curves (see Figure 31). 

3. Limiter and Nonlinear VCO 

Figure 35 gives the transient responses computed when 
both the nonlinear VCO model and the saturation effect were in 
the simulation program. Note that (1) the saturation portion of 
the curve is no longer a straight line, (2) the form of the re- 
sponse for (a) the linear AFC model, (b) the limiter, (c) the 
nonlinear VCO, or (d) the nonlinear VCO and the limiter are 

basically the same (Figures 31, 32, 34, and 35) and (3) the 

steady state frequency error for a given shift in ( v m ) re f is 
essentially the same for all the simulation models. The "final" 
frequency errors arrived at by the simulation model containing 
the nonlinear VCO and the limiter are tabulated for nominal loop 

gains ranging from 10 to 500 in Table IV. 


TABLE IV 

STEADY STATE FREQUENCY ERRORS WITH 
NONLINEAR VCO AND LIMITER IN THE AFC SYSTEM 


Ga in 

(Af ) 
o ss 

(k> 

(MHz) 

10 

.128 

25 

.054 

50 

.029 

100 

.014 

250 

.006 

500 

.003 
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GAIN, K, IS EVALUATED AT OPERATING 
POINT ON VCO CHARACTERISTIC 
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Figure 35 • AFC simulation transient response with nonlinear VCO 
and with ± 12 volt limiter. 



The table lists the steady state frequency error as a function 
of AFC loop gain for an open loop frequency shift of 0.8 MHz. 

In each case k is measured at the nominal VCO operating point. 

G. An Alternate Sampling System 

The feedback signal sampling circuit discussed thus far 
can be replaced by a simpler one with only a small loss in versa- 
tility. Deletion of the hold device in front of the RC charging 
circuit results in an AFC system that is somewhat slower in 
response at the higher gains, but is simpler to implement and is 
stable for all loop gains. The simplified system is not capable 
of "deadbeat" response, however. Consider the sampling scheme 
depicted in Figure 36 . 

1. Analysis 

The differential equation describing the linearized 
AFC system with the simplified sampler will be derived using 
Figures 26 and 36. An algebraic expression relating the sampler 
input, V^, to the sample -and-ho Id output, V^, and the system 
parameters derived from Figure 26 is 


V 3 = -kV 4 + D , 


(45) 


where the constants are 


k = kkk and 
A D1 v 


D = V + HkJf ■ nf - f + k (V - V )] 
bias A Dl L C R D v m c 


- w • 


(46) 




Figure 36. Alternate Sampler Model. 


Combining this result with a description of the sampler, the 
closed-loop AFC system differential equation can be written. 
The simplified sampler (Figure 36) is a charging RC circuit 
whose response is governed by 

de (t) 

(RC)~ dt + e o( t ) = e 2.( t ) ( 47 ) 

where e^(t) and e o (t) are the (sampled) input and output volt- 
age, respectively. Thus 

(RC) V 4 + V 4 = V 3 . (48) 

From consideration of the upper portion of Figure 26, 


V. = (f - f )/K + V - V , 

4 c o v me 


so that 


V. = -f /K 
4 o v 


(49) 


(50) 


Similarly, from the lower portion of Figure 26 one can deduce 
that 


V. = k * k m f + V, . + k [k (-nf - f ) - k (f ■ - f )] . 

3 A D1 o bias A L Dl R D D2 R o J 

(51) 
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Substituting Equations (49), (50), and (51) into Equation (25), 

and recognizing that V is constant during the sampling inter- 

m 

val, the differential equation governing the AFC system response 
during the sampling interval is 


r f c 

(RC) f + (k + 1) f = k iv-V +r £ -V, . 
o o vlmck bi 

v 

+ k A [ k I>l<" £ R + V + k D2 (f R ' V]} • 


bias 


(52) 


The system is seen to be stable for all positive values of gain, 
k, with the sampling-interval time constant being inversely pro- 
portional to the open loop gain for large gains, i.e., T = RC/(k+l). 

o -L 

The step response is always underdamped. Assuming that the sampling 

duration is t and that the sampling period is given by T, the 
° * 

effective time constant is given by 


t = RC/(k + 1) . (53) 

o 

For the typical parameter values used with the earlier sampler, 
the difference between the effective time constants for the two 
sampling schemes is negligible for loop gains below 500; that is, 
only when the deadbeat gain is approached does the earlier system 
show any advantage. Table V gives a comparison of the AFC system 
effective time constants for the two sampling schemes. The tran- 
sient responses of the two systems are essentially identical for 
normal AFC loop gains. The steady state responses . (sensitivities) 
of the two are also identical [see Equation (41)]. 


2. Simulation 


A simulation using the alternate sampler AFC system 


The "effective time constant" has significance when the loop 
gain is low enough that "several" sampling intervals are re- 
quired for recovery from a disturbance. 
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TABLE V 


COMPARISON OF AFC EFFECTIVE TIME CONSTANTS OF 

TWO SAMPLING SCHEMES FOR R = 20,000 OHMS 

C = 1 x 10" 7 FARADS, and t = 1.5 x 10 -6 SECONDS 

o 



T 


T 


Alternate 

Percent 

Flat-Top 

Gain 

Sampler 

Difference 

Pulse Sampler 

00 




-i 

00 


00 

-0.5 

2666.6667 

i — t 

O 

• 

1 

2667.1667 

0 

1333.3333 

0 

1333.3333 

1 

666 .6666 
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3.950 
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84.897 
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1.0658 

190.220 

.3598 
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along with the nonlinear VCO and limiter models was constructed. 

A listing of the simulation is included along with one simulation 
run (k = 250) in Appendix G. The response is not significantly 
different from that obtained with the flat-top pulse sample scheme 
(see Figure 35) for the lower gains (k = 10, 25, 50, 100). 

Figure 37 shows the response for AFC loop gains of k = 500, 750, 
1000, and 2000; for these gains the alternate sampler system is 
somewhat more sugglish than the sample-and-hold system when the 
error signal is below the limiting value. 


H. Conclusions 


A sample-data AFC system for direct FM transmitters has been 
modeled, analyzed, and simulated to establish the relationships 
between AFC system performance and the parameters of subsystems 
such as the VCO, the reference oscillator, the discriminator, and 
the sampling system. The analytic results obtained with the 
linearized closed loop difference equation include the following: 

(1) the range of loop gain, k, for which the system is 
stable and the form of the response as a function of 
the sampler characteristic parameter, C- , (Table I 
and Figure 28), 

(2) the time constant of the continuous envelope of the 
(discrete) transient response as a function of loop 
gain, k, and sampler parameter, C- [Equation (34) 
and Figure 29] , 

(3) frequency sensitivity; that is, the steady state fre- 
quency error as a function of shifts in the VCO center 
frequency, the input signal reference voltage, the 
reference oscillator frequency, and the amplitude 
dependence in the discriminator [Equation (41)], 

(4) a comparison of two sample-and-hold schemes. 

The linear version of the simulation was verified using the 
analytic results and was used to establish transient responses for 
various values of loop gain. Nonlinear models were then introduced 
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Figure 37* Alternate sampler AFC simulation transient response 
with nonlinear VCO and with ± 12 volt limiter. 




into the appropriate program blocks to simulate signal saturation 
in the feedback loop and to model the VCO characteristics more 
accurately. 

The nonlinear simulations indicate that the limiter does not 
destabilize the AFC loop, but only slows down the transient response 
until the error becomes small. This result is compatible with 
that obtained with the method of ’’describing functions” used in 
nonlinear control system analysis [8]. The latter treats a limiter 
as a variable gain which decreases as necessary to keep the out- 
put magnitude from exceeding a limiting value. The linear analysis 
shows that the AFC loop remains stable for decreasing positive 
values of loop gain. 

The simulation with a nonlinear VCO input-output characteristic 
similarly yielded a stable system. Action of the nonlinear VCO 
can be visualized as a variable gain. Note, however, that the 
apparent VCO gain increases for large negative perturbations in 
the VCO input voltage. Consequently, the system stability for 
high nominal loop gains is marginal. The simulation run made with 
both nonlinearities (VCO and limiter) included is also stable and 
well behaved. 

The two sample-and-hold schemes for feedback of the frequency 
error term yielded essentially equivalent performance for all except 
the highest AFC loop gains. Although the alternate scheme (instan- 
taneous error term driving the RC network) is theoretically more 
sluggish, the difference in response is negligible until the gain 
approaches that for deadbeat operation. The existence of a dead- 
beat gain seems to be the only advantage of the first scheme. The 
alternate scheme is much simpler to implement. 

Even with a limiter in the feedback path, the use of a relatively 
large loop gain seems desirable from a signal (as opposed to noise) 
point of view. The system still has high effective gain for low 
frequency errors producing small steady-state frequency errors. 

The transient response for small perturbations of parameter values 


91 



is also much faster with high gains. From a noise point of view, 
on the other hand, increased gain is detrimental. Noise effects 
should be studied using simulation and/or analytic results before 
raising the loop gain of any existing AFC systems to obtain the 
faster response and smaller steady state errors. 
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IV. DETAILED CIRCUIT SIMULATION: STIFF SYSTEMS 


A. Introduction 

The digital simulation of a linear dynamic system is often 
complicated by the existence of a large spread in the eigenvalues 
of the system. Since the step size in the numerical integration 
needs to be small with respect to the shortest time constant of 
the system while the range of time over which integration is 
desired is usually greater than the largest time constant of the 
system, a large spread in eigenvalues leads to an excessive number 
of steps in the numerical integration. Simulations of such stiff 
systems suffer from large round-off errors and excessive computer 
time requirements. In those cases where the transient response 
associated with the short time constant modes is of no interest 
(having decayed to an insignificant level almost immediately) , 
one would like to "remove” those fast modes from the system model. 
This section describes the development and application of a method 
for removing such roots from the equations comprising the model. 
This permits approximating a linear dynamic system by a set of 
differential equations of lower order not containing the fast 
eigenvalue(s) . 

Following the development of the theory of "root removal," 
a simulation program was developed which uses the root removal 
technique to efficiently simulate linear stiff systems. The 
simulation, known as the "Approximate Linear System Analysis 
Program" (ALSAP) , accepts the system differential equations in 
state variable form, approximates them by a lower order system 
of equations without the fast eigenvalues, and predicts the sys- 
tem response by numerically integrating the approximate system 
equations. The construction of ALSAP is described, along with 
its application to a thirteenth order system representing an 
RLC bandpass filter. The system response predicted by ALSAP 
compared favorably with the response of the actual filter as 
obtained in laboratory measurements. 
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B. Derivation of Root Removal Technique 

Assume that a linear, time-invariant system is described in 
the usual state variable form, 

x = Ax + Bu . (54) 

One approach for "removing" a particular mode would be to first 
transform the system so that the matrix A is in Jordon canonical 
form. The scalar differential equation corresponding to the fast 
mode could then be assumed to have an approximate solution equiv- 
alent to an instantaneous response. The remainder of the scalar 
differential equations would not contain the fast root and could 
be numerically integrated if desired. The solution vector would 
then be transformed back from the canonical frame to the original 
coordinate frame. Note that the matrix required for the transfor- 
mation has all of the system eigenvectors as its columns and that 
the inverse of this matrix is required [7]. 

The above procedure is not practical. However, a related 
procedure has been developed which requires only the fast eigen- 
value and its single associated eigenvector in order to remove a 
fast root. The procedure does not require matrix inversion. The 
resulting approximation is analogous to replacing the transfer 
function [l/(s+a)j by the gain [l/a] in a system for which the 
time constant (-1/a) is much faster than any of the other time con- 
stants of the system or components of the system input signal. 

The resulting approximate system is related to that which Chidambara 

[9] developed for control system synthesis. The digital simulation 
application is not sensitive to the points of criticism of Davison 

[ 10 ] . 

Let the matrices in Equation (54) be partitioned to isolate 
the last element of the state vector: 
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n-1 


n 


L 11 


l> 


12 


nn 


n-1 


u . (55) 


It is desired to find a simple n-by-n transformation matrix, P, such 
that the state equation [Equations (54) or (55)] expressed in a new 
coordinate system has the fast root isolated for approximation. 

Let the transformed state vector, 2> be defined by 


x = Py 


(56) 


The resulting (transformed) state equation is given by 


2 = (P“ l AP) 2 + (i’” 1 B) u . 


(57) 


Let the fast root (eigenvalue) to be removed be \, and let its 

eigenvector (n-by-1 matrix) be v. Assume that the eigenvector 

th 

has been normalized such that its n component is unity. A satis- 
factory transformation matrix, P, will be shown to be formed by 
replacing the last column of an n-by-n identity matrix by the eigen- 
vector, v. Partition the vector into two parts such that its first 
(n-1) components form the vector j> . That is, such that 


* ■ [ 4 -] • 


(58) 


Then the transformation matrix can be written in partitioned form as 

£ 


P = 


_I n-l ! 


o' 


(59) 
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T 

where 1^ ^ is an (n-1) by (n-1) identity matrix and 0 is a 1 by 
(n-1) matrix of zeros. Note that one of the features of this partic- 
ular transformation matrix is that its inverse can be formed by 
simply reversing the signs of the top (n-1) elements in the n'^ 1 
column of the matrix P; that is, 



The term (P ^AP) which appears in the transformed state equation 
[Equation (57)] can best be evaluated by multiplication of the par- 
titioned matrices: 




^9 1 I ^01 £ + a _„ 


It is desired to show that the upper-right-hand submatrix of P ^AP 
is zero and that the lower-right-hand submatrix is the fast eigen- 
value, To show this, recall the definition of eigenvalue and 

eigenvector: 


Av = \v . 


( 62 ) 
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Writing this definition in partitioned form, one has 


f A n : 

_^12 



= X 

£ " 

A. , i 

a 


1 


21 * 

nn 



1 


which can be separated into two equations as follows: 


( 63 ) 


and 


A 11 £ + A 12 = 

A 01 £ + a = X . 

21 nn 


(64) 

(65) 


Substitution for X from Equation (65) into Equation (64) gives (since 
\£ = £\) 


A 11 £ + A 12 ■ £ (A 21 £ + V • 

A u £ + A 12 - £ A 21 £ - £ a n n * 2 • 


( 66 ) 

(67) 


Thus the 
Equation 
the fast 


upper right submatrix of P ^AP is indeed zero. Furthermore, 
(65) shows that the lower right element of P ^AP is equal to 
eigenvalue, Thus 


P _1 AP = 


A, „ - £ A_ , 0 


11 


21 I - 


21 


( 68 ) 


Note that this matrix is easily constructed from the matrix A, the 
fast eigenvalue X, and the upper portion of the associated eigenvector; 
note further that matrix inversion is not required. 

The matrix (P ^B) is also required in Equation (57): 



( 69 ) 
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Before substituting the expressions for P ^AP and P into Equation (57) 
partition the new state vector, £, into two sub-vectors with the second 
being a scalar. That is. 





(70) 


Substitution of Equations (68), (69), and (70) into Equation (57) allows 
the transformed state equation to be written in the following form: 


• 

r 


" A ll - £ A 2l 

' “1 


r 

-L 

f B l - £B 2l 

y n 

L J 


I— 1 

CM 

< 

1 

1 \ 
I 


.v 

T 

1 

CM 

PQ 

J 


Note that the variable, y , does not enter into the r equation; the upper 
(n-1) scalar differential equations can be solved independently of the y^ 
equation. It can be shown that the fast root isolated in the y^ equation 
does not appear in the r equation. 

The transformed state equation can be separated into two state 
equations : 


- = ( A ll ■ £ A 21^ - + ^ B 1 " £ B 2^ - 

and 

y n = A 21 £ + Xy n + B 2 £ ' (73) 

Equation (72) is a state equation of order (n-1) whose eigenvalues are 
those of the original system [Equation (54)] except that the fast eigen- 
value, \, is absent. Thus Equation (72) can be numerically integrated 
with reasonable simulation time requirements. To reconstruct the original 
state vector, the value of y (t) as well as .r(t) is required. Keeping 
in mind that jr(t) can be evaluated independently of y^(t), rewrite 
Equation (73) with r(t) considered as an additional input. 

y n = X y n + (A 21 r + B 2 u). (74) 
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While this differential equation could be solved explicitly, it is 
useful to approximate its solution by assuming that |\| is suffi- 
ciently large (fast) that the response in y^(t) to a change in r(t) 
or u ( t) is essentially instantaneous. This is equivalent to replacing 
the transfer function [l/(s+a)] by the gain [l/a]; this is often done 
in frequency domain studies when the root is much faster than input 
variations or other roots of the system. The resulting approximation 
is 

y n «"ir (a 2 i — + b 2 — } * (75) 


C. Reconstructing the Original State Vector 

Assuming that the _r equation has been successfully integrated 
and that the y^ equation has been approximated as above, the original 
state vector, x> can be reconstructed as 


x = Py = 


^-i! 

£ 


£ 


r + py 
— ^n 

LV : 

1 


_ y n _ 


y n 


(76) 


Substitution of Equation (75) into Equation (76) gives 


h-l ‘ £ A 21 /X 


” £ B 2 /\ 

_ " a 2 ]/ x 

r - 

b 2 /x 


(77) 


The initial conditions for the components of jc can be expressed in 
terms of the initial conditions of the components of x as follows: 

r.(o) = x. (o) - p. x (o) . (78) 

l l in 
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Recapitulating, for a set of state equations as in Equation (54) 
containing one fast root, a lower order system [Equation (72)] which 
does not contain the fast root can be numerically integrated using 
the initial conditions as given by Equation (78); the original state 
vector, x(t) > can then be determined by applying Equation (77)'. 

These steps are well-suited for inclusion in a simulation and the 
computer time required for the transformation will be small in com- 
parison with the savings in integration time. 


Example : 


,rd 


Consider the following 3 order system with scalar control: 


x. 




— — 


- — 


— -i 

-3 9 8 


X 1 


0 

1 -11 -8 


X 2 

+ 

0 

-1 -9 -12 


x~ 


1 



3 


L J 


u. 


(79) 


The roots of the system are -2, -4, and -20; assume that the "fast* 

root, \ = -20, is to be removed. The corresponding eigenvector is 

T 

found to be v = (-1, 1, 1) so that 


£ = 


-1 

1 


(80) 


The necessary terms in Equation (72) are evaluated as follows: 






r_ _ 

r 


r- — 


-3 9 


-l 

L' 1 " 9 . 


O 

I 

<A U - £ A 21 > > 

_ ! -11 

- 

. i. 


A 

_ 2 - 2 J 


(81) 


(B, - £ B-) - 


1 

_ 0 , 


’ -l 


l" 


- 


II 

i — l 
rH 
l 1 


o 7 

1 


l 


-1_ 


(82) 


Repeated application is required when more than one fast mode is 
to be removed. 
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The reduced-order system then is 



The initial conditions are found from Equation (78) to be 


(83) 


r^(o) = x 1 (o) + x 3 (o) 

and 

r 2 (o) = x 2 (o) - * 3 (o) , (84) 


and the original state variables are expressed in terms of the re- 
duced system variables by applying Equation (77) . 


’ X 1 


T .05 

.45* 


r l 


~-.05~ 

X 2 

= 

-.05 

.55 


r 2 

+ 

.05 

X 3 
L J 


-.05 

-.45 




.05 


Thus the solution to Equation (79) can be approximated by simulating 
the reduced system described by Equations (83) and (84) and substitut- 
ing the results into Equation (85). This procedure is straightfor- 
ward to implement in a digital simulation. 

D. Practical Application: RLC Bandpass Filter 

The CCS telemetry bandpass filter which was analyzed earlier in 
this report from a frequency domain point of view and whose schematic 
is given in Figures 1 and 2 has been simulated using ALSAP. The state 
variables were chosen as the capacitor voltages and inductor currents, 
except that the (excess) current through the mutual inductance was 
not included. The specific relationships between the state variables 
and the circuit voltages and currents are specified in Table VI. 
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TABLE VI 


PHYSICAL INTERPRETATION OF THE BANDPASS FILTER 
STATE VARIABLES 


State Variable 

Physical Interpretation 
(volts, amps) 

x i 

Voltage across 

c i 

X 2 

ii it 

C 2 

X 3 

ii ii 

C 3 

X 4 

ii it 

C 4 

X 5 

ti ii 

s 

X 6 

ii ii 

C 6 

X 7 

it it 

C 7 

X 8 

Current through 

l ti 

X 9 

ti ti 


X 10 

it it 

h 

X 11 

ii n 

L 3 

X 12 

it it 

L 4 

X 13 

ti m 

L 5 


Note: The circuit elements referred to here are shown in Figures 1 

and 2, Positive currents are assumed flowing to the right 
or down and voltages are considered positive at the left or 
at the top of the element. 

Also, the component nomenclature for the ALSAP simulation is 
related to the component nomenclature shown in Figure 2 by 
the following: LT1 = Lll, LT2 = L13, LM = L12 , RT2 = R13,. 

RL2 = R4, RL3 = R5 , RL4 = R6, RL5 = R7 , and the double sub- 
scripted C values shown in Figure 2 were not used in the 
ALSAP simulation 






The state equation was written using Bashkow's A matrix method [ 11] 
and is of the form x = Ax + Bu, where x is the thirteenth order 
state vector, u is the input voltage, A is the thirteen-by- thirteen 
state distribution matrix and B is the thirteen-by-one control dis- 
tribution matrix. The resulting A and B matrices are specified 
in terms of the circuit elements in Table VII. The thirteen eigen- 
values of the matrix A are given in Table VIII. 

The nominal operating frequency of the filter is 1.024 MHz. 

It can be seen from the system eigenvalues in Table VIII that two 
of the eigenvalues, \ ^ and , are trivially fast with respect 

to the other eigenvalues. Using AI5AP, these two fast roots were 
removed; the simplified eleventh order system was numerically 
integrated; and the original state variables were reconstructed 
for printout. While numerical integration of the original state 
equations would have required a prohibitive amount of machine time, 
the time required for ALSAP to transform and integrate the system 
equations for a reasonable response time was under one minute. 

Figure 38 shows the filter response to a one MHz input as deter- 
mined by ALSAP. Laboratory tests on the actual filter confirmed 
the simulation results, showing that the root removal approximations 
were valid. Appendix H gives a listing of ALSAP. 

E. Complex Pairs of Fast Roots 

When the fast root to be removed is complex, the same procedure 
as used for removal of a real root can be used provided the simula- 
tion uses complex arithmetic. However, complex roots always appear 
in complex conjugate pairs and any complex roots that are "fast' 1 
must be removed in pairs. The necessary equations have been derived 
for removing a fast complex pair of roots in one step without using 
complex arithmetic in the simulation. These equations, which are 
somehwhat more complicated than those for single root removal, are 
given in Appendix I. 
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TABLE VII 


NON-ZERO ELEMENTS OF STATE AND CONTROL DISTRIBUTION MATRICES 
FOR BANDPASS FILTER SHOWN IN FIGURES 1 and 2 


All elements of the A and B matrices are zero except for the follow- 
ing elements: 

A( 1 , 10) = 1/Cl 
A(l, 11) = 1/Cl 
A(2 , 10) = 1/C2 
A(3, 11) - 1/C3 
A (4, 11) = 1/C4 
A (4, 13) = -1/C4 
A(5 , 11) = 1/C5 
A(5 , 12) = -1/C5 
A (5 , 13) = -1/C5 
A(6,6) = -1/ (R3 *C6) 

A(6 , 13) = 1/C6 

A(7 , 7) = -1/ (RL*C7) 

A(7, 13) = 1/C7 

K1 = (LT2 + LM)/(LM-[LT1-LT2 + LM° (LT1 + LT2 ) ] ) 

K2 = 1/(LT1-LT2 + LM*[LT1 + LT2]) 

A(8 , 8) = -LM-R1-K1 

A(8,9) = - (RT2 + R2) • (LT1 + LM) *K1 + (R2 + RT2)/LM 

A(8, 10) = R2 • (LTl + LM) »Kl - R2/LM 

A(8, 11) = R2 • ( LT 1 + LM) *K1 - R2/LM 

A (9 , 8) = -m*Rl'K2 

A(9,9) = - (RT2 + R2) • (LTl + LM) -K2 

A(9 , 10) = R2 • ( LT 1 + LM) *K2 

A(9, 11) = R2 • (LTl + IM)*K2 

A(10, 1) = -1/L2 

A(10,2) = -1/L2 

A(10 , 9) = R2/L2 

A (10, 10)= - (R2 + RL2)/L2 

(CONTINUED) 
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TABLE VII (CONTINUED) 


NON-ZERO ELEMENTS OF STATE AND CONTROL DISTRIBUTION MATRICES 
FOR BANDPASS FILTER SHOWN IN FIGURES 1 and 2 


A(10, 11) = -R2/L2 

A(ll , 1) = -1/L3 

A(ll,3) = -1/L3 

A(ll,4) = -1/L3 

A(ll,5) = -1/L3 

A(ll, 9) = R2/L3 

A(ll, 10) = -R2/L3 

A(ll , 11) = - (R2 + RL3)/L3 

A (12 ,5) = 1/L4 

A(12 , 12) = -RL4/L4 

A(13,4) = 1/L5 

A(13,5) = 1/L5 

A(13, 6) = -1/L5 

A(13,7) = -1/L5 


A(13 , 13) = -RL5/L5 


B (8) = LM-K1 
B (9) = LM-K2 




TABLE VIII 


EIGENVALUES OF THE BANDPASS FILTER A-MATRIX 


X 2 = (-6.11 X I0 4 ) + j(O) 

\ 3 = (-4.18 X 10 5 ) + j (-6.31 X 10 6 ) 

\ 4 = (-4.18 X 10 5 ) + j (6.31 X 10 6 ) 

\ 5 = (-1.87 X 10 5 ) + j (-5.84 X 10 6 ) 

\ 6 = (-1.87 X 10 5 ) + j (5 .84 X 10 6 ) 

X ? = (-7.81 X 10 5 ) + j(O) 

X g = (-2.44 X 10 5 ) + j (-6 . 74 X 10 6 ) 

X 9 - (-2.44 X 10 5 ) + j (6.74 X 10 6 ) 

X 1Q = (-1.03 X 10 5 ) + j (-1.77 X 10 7 ) 

X n = (-1.03 X 10 5 ) + j (1.77 X 10 7 ) 

X 12 = (-2.49 X 10 8 ) + j (0) 


X B = (-1.96 X 10 9 ) + j(0) 
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F 


Conclusions 


A method has been developed for efficient digital simulation 
of stiff linear dynamic systems by the "removal" of insignificant 
eigenvalues. The method has been implemented in the digital simula- 
tion ALSAP and has been successfully applied to the simulation of 
both passive and active networks. The computer time required for 
the simulation of a typical linear network using AliSAP is signifi- 
cantly less than that required for available general circuit analy- 
sis programs. The general programs, however, are capable of gene- 
rating internally the governing differential equations from a 
topological description of the network while ALSAP currently re- 
quires that the engineer write the network state equations. The 
inclusion of an option allowing the use of the ALSAP technique 
with stiff linear systems in a general computer aided design pro- 
gram would result in a more useful simulation facility. 
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APPENDIX A. 

Subroutine FOFT - for evaluation of 
the complex time function 
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APPENDIX B 


Subroutine TLMFLT - the evaluation of the CCS telemetry 
bandpass filter transfer function. 


113 



GUPROU TlNF TLVFLfU.N.NtsM l 

C THIS S'JRROUT i Vt CALCiJLATrs THE TRANSFER FUNCTION l/F urtE CCS TELEMETRY 

C tsANuPAJS FlLTlk* NORMAL I 2 1. G IT TO UN ITT AT IT* PtA< V ALUF, AND APPLIES 

C IT TO THE INPUT fkr'.ULNCtf A- 4 Ay A. 

0 t» O 9 0 1> Q fltt-tt* 0 fr » * 0 » a C <>■{>{} 

COMPLEX 1 1 1 i*12*Zn./2.?'5.£* # /3,26**7,4b.Z9.21u«2M*iA2.2a.ZC*ZDD. 

1ZL./F ,ZN*7UtZf>l •7.0?,ZL»3 , 2 Da »z 1 5 ,206 *ZD 7 ,2D8 *209 ,H , A * e 3 6 ) f 3 , 256 ) 

RtAL LI ) »Ll?9LlitL2.L2*LA*L5 
N2 = N / 2 
NN = N - 1 
T 2 = 0. 

C THE FOLLOWING INSTRUCTIONS SET THE FILTER ELEMENT vALOES 

R2 = 5.1E3 
R3 * 4 • 7E 3 
RL = 6 • 2E2 
Cl » 4.80F-10 
C2 = 1.368E-9 
C3 = 5.10E-10 
C4 = 2.2E-9 
C5 = 3.42E-10 
C6 = 5.7E-10 
C 7 = 1.8E-9 
C22 = 2.07E-11 
C33 * 1.65E-11 
C 5b = 1.65E-11 
Lll = -5 • 358E-6 
L 12 = 6 • 06 E-6 
L 1 3 = b • 19AE-5 
L2 = 1.08E-5 
L3 = 5.1 E-5 
LA = 1.6E-5 
L5 = 5.1 E-5 

C DO LOOP 10 CALCULATES THE FILTER TRANSFER FUNCTION 

DO 10 I = 2*N 
11=1-1 

FI = 1. / (N3IT * l.E-6) 

F = FI * I I 
W = 6.2831853 # F 
WSQ = w ** 2 

211= CMPLX ( 5 1 • . W * Lll ) 

212 = CMPLX (0. * W * L12 > 

213 = CMPLX(7o78* W * L13) 

22 = CMPLX ( R2 9 0* I / CMPLX<1.» W * R2 * C22) 

23 = CMPLX ( 0. 9-1. / IW * Cl)) 

24 s CMPLX (1.70. IW * L2 - 1. / CW * C2) M 

25 = CMPLX ( 6 .0 1 • W * L3 ) / CMPLXMl. - WSQ * L3 * C331 * 

1(W * 8.01 * C33M + CMPLX (Oa* - l./tw * C3)) 

26 = CMPLX(2o51 * W * L4 ) / CMPLXMl. - WSQ * L4 # C5l. 

1 (2.51 * W * C5>) + CMPLX ( 0 • * -1. / (w * C4>) 

27 = CMPLX ( 0 o 0 1 t W * L 5 1 / CMPLXMl*- wsG * L5 * C55»* 

1 ( W » 8 . U 1 * C 5 5 M 

28 = CMPLXU.E7, O.I 

29 = CMPLX ( R 3 » 0.) / CMPLXtl.* W * R3 * C6) 

210 = CMPLX(RL*0.) / CMPLXI1.* W * RL * C7) 

2 A 1 = 211 + 212 

2A2 = Z 1 2 + 213 + 22 

ZB = 22 + 23 a Zh 

ZC = 24 + 25 + 26 

ZD = 26 + 27 + 20 

2E = Z8 + 29 + 210 

2F = 2 A 1 * 2A2 - (212 ** 2oI 

2N = 212 «■ 22 » 24 * 26 » 20 * 2lu 

2D 1 = 2D * 2E * (24 ** 2.) 

2D2 = 2B * 2E * (26 <>* 2d 

2D3 = (24 2 o ) * (28 *© 2. ) 

ZD4 = 23 # 2C * 2D * ZE 
2D5 = ZB « ZC ^ (28 ** 2.) 

2D6 = ZA1 * C Z 2 ** 2. I 
2D7 = 2C * 2D * ZE 

ZDS = ZC o ( 28 ** 2* ) 

209 = ZE * ( Z 6 ** 2* J 

2DD - If o (-ZD1 - 202 «■ 203 + 2 D4 - 705) 

1 ♦ 2D6 * ( -ZD7 ♦ ZDS + 209) 

H = 2N / 2 DD 
T 1 = CA3S(H> 

IF( II 0LT0 T 2 ) GO TC 15 
T 2 = T 1 
15 CONTINUE 
10 B ( I ) = H 

WR I TE ( 6 * 1002 ) 

1002 FORMAT ( 2 HI > 

WRITE (6,1001) ((1*8(1)), I = 2.N) 

1001 FORMAT ( 1 H . 1 6« 2E20.6 ) 

COR = I, / T2 

C DC LOOP 25 NORMALIZES THE PEA* vAlUE uF THE 1 RAnSFEk FUnCTIu* TO 1*0 AnD 

C APPLIES IT TO THE INPUJ ARRAr A 

DO 25 I = 1 , N 
6(1) = COR o B ( I) 

25 A ( I ) = A ( J ) o (HI) 

WRITE ( 69 1020) 

1020 FORMAT ( 1 HO 9 2 X • 20HNORMAL 1 2 A T I ON FACTOR) 

WKITr ( 691015 ) COR 
1015 FORMAT ( 1 H+ 9 ? 3 x 9 L . T 0 . 6 ) 

Ml URN 

p % • 
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APPENDIX C. 

Main program S-band calculations. 
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COMPLEX AC 819 2 ) 

DIMENSION I BUF ( lOOOO ) 

CALL PLOTS ( I6UF ( 1 ) ,10000,2 ) 


I GAM = 13 
BETA = 1.0 
FC = 1 © 0E6 
FMOD = 1.36533 33E4 
F PCM = 4.9834666E5 
F SC = 1.024E6 
PI = 3.1415927 
P 1 2 = 6.2831853 
NBIT = 146 

N = 2 ** I GAM 

W = PI 2 * FSC 

WM = PI2 * FMOD 

WPCM = PI2 * FPCM 

DELT = (2. / ( FMOD * N ) ) 

DELF = FMOD / 2. 

DO 1 I = It No IT 
C ( I ) = 1 

1 CONTINUE 

DO 2 I = 2 * Ncs I T » 2 

2 C ( I ) = -1 

CALL TIMFCN(A,C.N,DELT , W * W M , vi P C M ) 
CALL FOFT (A, IGA M, BETA) 

CALL F F T ( A 9 I G A M s - 1 ) 

CALL L F 0 L D ( A 9 N ) 

CALL FFTPLT (A, N, DELF ) 

CALL PLOT ( 0 • 9 0 • *999 ) 

STOP 

END 
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APPENDIX D 


Subroutine TIMPLT - Calcomp plot routine for 
plotting baseband time functions. 



SUBROUTINE T!MPLT(A,N»PRD) 

c *«**•#*»•**»*•*«*»•##»*»*»»»«»•«*•••••« 

COMPLEX A < 1 > 
call FACT OR (0*4) 
call PLOT ( 0. .-20. ,3 ) 

CALL PLOTt 12. .0.9-3) 

CALL PLOTf-lO. ,-14.,3 ) 

(;«■*«■****#*#****■#**###*#*#**##*#»-».####### 

C DO LOOP 3 DRAWS THE uORDER OF THE PLOTTING AREA 

C**»»**«*«*-**«##*»*4»***»****»*#ft»»««««« 

DO 3 I - 1*2 

CALL PLOT( -10. ,0.»2 ) 

CALL PLOT! 10. .0. ,2 ) 

CALL PLOTf 10. *-14. *2 ) 

CALL PLCTI-10.01.-1A. ,2 I 
CALL PLOT( -10.01*0.01 *2 ) 

CALL PLOTUO. 01 *0.01.2) 

CALL PLOT ( 10.01 , -1 4 • Cl «2 ) 

3 CALL PLOT(-10. ,-14.01,2 > 

(!*******#*•****»***•#»*»•*•»•«**•«• 

C DO LOOP 10 TICKS THE POSITIVE HALF OF TH£ LEFT HAND rORDER 

C***« *****»»*»»***.*#•»*****.. *»*«**»•«• 

DO 10 I * 0,16 

Y = - 7. + 0.4375 * I 

IF < MOD ( 1,5) .EQ. 0 1 GO TO 6 
CALL PLOT ( -10. 1 ,Y*3 I 
GO TO 8 

6 CALL PLOT( -10.2 ,Y,3) 

8 CALL PLOT! -10. ,Y,2 1 

10 CONTINUE 

C**#*a««e*»*««»0**#*»****»«*»»*»ft»»»*»«» 

C DO LOOP 11 TICKS THE NEGATIVE HALF OF THE LEFT HAND FORDER 

C-fro***#*****#***-*#*#******-*#.**#*#****#*# 

DO 11 I = 1*15 

Y * - 7. - 0.4375 * I 

I F ( MOD ( I ,5 1 .EQ. 0) GO TO 7 
CALL PLOT < -10. 1 . Y * 3 ) 

GO TO 9 

7 CALL PLOT ( -10.2 *Y*3 ) 

9 CAlL PLOT(-10. ,Y,2 ) 

11 CONTINUE 

C*#«#******************#««***»»»»*#»«4»» 

C DO LOOPS 12 AND 13 AND THE FOLLOWING INSTRUCTIONS PROVIDE THE CALIBRATION 

C AND LABELING OF THE LEFT HAND BORDER 

C****#*****»#*-&***##*###*##**. < ,**4fl.»^{>000 

DO 12 I = 0*3 

Y = -7.105 + I * 2.108 

12 CALL NUMBER(-10.88,Y, .21.1.0*1 *0.*1) 

DO 13 I = 1,3 

Y = -7.105 - I * 2.188 

13 CALL NUMBER(-11.09*Y, .21.-1. 0*1. 0.*i) 

CALL SYM30L (—11 .03,-3 .5, *21 »17hamPL! T uDE < VOLTS ) .90. , 1 7 I 
C**-*****#**-*.**-**}-* »**#**•«■ «**»«. «»»•*»**«• 
C DO LOOP 25 PROVIDES THE TICK MARKS ALONG THE BOTTOM FORDFR 

C************#*********.***#**** • ****«•* 

DO 25 1 = 0,100 
X = -10. + 0.2 * I 
I F { MOD (1.5) .EQ. 0) GO TO 26 
CALL PLOT (X ,-14.1,3 ) 

GO TO 28 

26 I F ( MOD ( 1*10) .EQ. 0) GO TO 27 
CALL PLOT (X. -14. 16. 3) 

GO TO 28 

27 CALL PLOT (X. -14.2*3) 

20 CALL PLOT ( X .-14. .2 ) 

25 CONTINUE 

C DO LOOP 20 PROVIDES THE TICK MARKS ALONG THE TO p BORDER 

C**«-*C***4*#***#*4.*4*4. # 4 # ** # #** ### h ### 

DO 20 I = 0,100 
X = -10o ♦ 0.2 * I 
IF (MOD (1.5) .EQ. 0) GO TO 21 
CALL PLOT (X ,0. 1 ,3 ) 

GO TO 22 

21 I F ( MOD ( I . 1 0 1 .EQ. 0) GO TO 23 
CALL PLOT (X ,0.16,3 ) 

GO TO 22 

23 CALL PLOT ( X ,0.2,3 ) 

22 CALL PLOT (X ,0.0,2 ) 

20 CONTINUE 
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nnnn nnn onn 


DO LOOP 35 TICKS THE POSITIVE HALF OF THE RIGHT HAND BORDER 
»•***««**»****##*»»••* + ««#*«#«#***»***» 

DO 35 I = 0,16 
Y * - 7. + I * 0.4375 
I F ( MOD (1*5) .EG. 0) GO TO 36 
CALL PLOT ( 1 0. 1 , Y ,3 ) 

60 TO 37 

36 CmLL PLOT ( 1 0. 2 * Y *3 ) 

37 CALL PLOT i 10. ,Y,2} 

35 CONTINUE 


DO LOOP 40 TICKS THE NEGATIVE HALF OF THE RIGHT HAND 30RDER 


DO 40 I* 1,15 

Y = - 7. - I * 0.4375 
I F ( MOD (1*5) .EO. 0) GO TO 41 
call PLOT (10,1, Y, 3) 

GO TO 42 

41 CALL PLOT ( 10. 2 »Y« 3 ) 

42 CALL PLOT ( 1 0. , Y ,2 ) 

40 CONTINUE 

THE FOLLOWING INSTRUCTION , DO LOOP 50. AND THE INSTRUCTION FOLLOWING 
DO LOOP 50 PROVIDE THE CALIBRATION OF THE LOWER BORDER. 

CALL NUMBER (-10. 105, -14. 5,. 2 1,0. ,0.,-l ) 

DO 50 I = 1,9 

X * -10.21 +1*2. 

50 CALL NUM8ER(X,-14.5*.2l ,20. *1 

CALL NUMBER (9.685»-l4.5t.21,200.*0**-l) 

C 

CALL SYMBOL {-2.*-14.9».21»19H TIME ( M I CROSE CONDS ) *0 . , 1 9 > 

C 

C CHANGE OF ORIGIN 

C 

CALL PLOT(-10. ,-7. ,-3 ) 

c i »«« ft *»»»««»»»*»** t »««»«««« t «»*««»« t t 

c plotting the data with a straight line between point* 

C**»**«ft«ft««»*»»««»»«ftft*««««ft«ftftftftft«ft# 

DO 60 I « 1 , N 
11=1-1 

X = { I I * PRD* (0.2E6 ) ) / N 
I F ( X .GE. 20. ) GO TO 71 

Y * 2.1875 * REAL! A ( I ) ) 

60 CALL PLOT(X,Y,2) 

DO 70 I * 1,N 
II = I - 1 

X = (II * PRD* I0.2E6)) / N + PRD * 0.2E6 
I F ( X .6T. 20. ) GO TO 71 

Y = 2.1875 * REAL I A ( I ) I 

70 CALL PLOT(X,Y,2) 

c ft ft » ft ft ft »« « »« « « t »« t « « « » * ft ft ft ft ft • t « « ft * « ft ft » 

C RETURNING THE PEN TO THE LOWER EDGE OF THE PAPER 

C « ft ft « • * ft « « « « « * * ft « ft « * « » * • » t « » * « ft » « « » « « « 

71 CALL PLOT! 26. *-13. f-3 > 

C 

WRITE (6* 1035) 

1035 FORMAT < 1H1 ,2X, 14HPL0T COMPLETED) 

RETURN 

END 

C * * * * * * * * * * * # * * * # t * * * * » t * « * t * « . t * 
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APPENDIX E 


FORTRAN program for linearized AFC system simulation. 
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FKcCJUENCY CONTROL SIMULATION 


LOOP GAIN 100. 


STuP 

r 

VM 

V3 

PO 


MCStC 

VOLTS 

VOLTS 

HZ 

1 

63. 

.0000 

-7.0000 

•225000409 

2 

12d. 

.0000 

-7.0000 

.225000409 

3 

169, 

.0000 

-7.0000 

,225000409 

4 

2b<i. 

.0000 

-7.0000 

.225000409 

5 

315. 

• 0000 

-7.0000 

,225000409 

6 

37d, 

.uooo 

-7.0000 

.225000409 

7 

441. 

.0000 

-7.0000 

.225000409 

8 

504. 

.U0O0 

-7,0000 

.225000409 

9 

5b7. 

.0000 

-7.0000 

•225000409 

10 

63u • 

1,0000 

93,0000 

,226400+09 

U 

693. 

1,0000 

85,5027 

.226295+09 

12 

7bo. 

1.0000 

/a. 5731 

.226198409 

13 

819, 

1.0000 

72,1684 

,226108+09 

14 

862 . 

1.0000 

66.248o 

.226025+09 

lb 

94b. 

1 . 0000 

60.7770 

.225949+09 

16 

lOUd, 

1.0000 

55,7199 

.225878+09 

17 

1071, 

1.0000 

51,045b 

.225813+09 

18 

1134, 

1.0000 

46.7251 

.225752+09 

19 

1197, 

1.0000 

42.7320 

.225696+09 

20 

1260, 

1.0000 

39.0411 

.225645+09 

21 

1323, 

1.0000 

35.6299 

•225597+09 

22 

1360 o 

1.0000 

32.4767 

.225553+09 

£.5 

144 V, 

1.0000 

29.5626 

.225512+09 

24 

1512. 

1.0000 

26.8689 

.226474+09 

25 

157b. 

1.0000 

24.3793 

.225439+09 

2b 

1636, 

1.0000 

22.0781 

.225407+09 

c 7 

17U1 o 

1.0000 

19,9513 

.225377+09 

*8 

1764, 

1.0000 

17,9854 

.225350+09 

29 

1827 , 

1.0000 

16,1664 

.225324*09 

30 

189u g 

1.0000 

14.4890 

.225301+09 

31 

1953. 

1.0000 

12.9369 

•225279+09 

32 

10 lo. 

1.0000 

11,5021 

,225259+09 

33 

2079, 

1.0000 

10.1761 

,225240+09 

34 

2142, 

1.0000 

8.9504 

.225223+09 

35 

220b , 

1.0000 

7.8177 

,225207*09 

36 

2266. 

1.0000 

6,7706 

.225193+09 

J 7 

2331, 

1 .0000 

5.8029 

.225179+09 

08 

239h 0 

1.0000 

4.9084 

,225167+09 

J 9 

2457 . 

l.OQOO 

4.0816 

.225155+09 

4U 

252o, 

l.uooo 

3.3174 

.225144+09 


UELF 

HZ 

0 * 

0. 

0 . 

0. 

. 0 . 

o* 

~ 0 . 

0. 

0. 

1400000. 

129S038. 

1198024. 

1108358. 

1025480. 

948878. 

878078. 

812638. 

752152. 

696248. 

644576. 

596813. 

552674. 

511876. 

474164. 

439310. 

407094. 

377318. 

349796. 

324358. 

300846. 

279116. 

259030. 

240466. 

223306. 

207448. 

192788. 

179240. 

166718. 

155142. 

144444. 
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FREQUENCY control simulation 




LOOP GAIN 100, 



STtP 

T 

Vrt 

V3 

FO 

UELF 


MCScC 

VOLTS 

VOLTS 

HZ 

HZ 

41 

2563, 

1.0000 

2.6111 

,225135+09 

134556. 

42 

264o, 

1.0000 

1,9584 

.225125+09 

125418. 

43 

27u9, 

1.0000 

1.3550 

•2251 17+09 

116970. 

44 

2772. 

1.0000 

,7973 

,225109+09 

109162. 

4b 

283b. 

1.0000 

,2819 

.225102+09 

101946. 

46 

289d. 

1.0000 

-.1946 

.225095+09 

95276. 

47 

29ol, 

1.0000 

-.6349 

.225089+09 

89112. 

46 

30^4, 

1.0000 

-1,0419 

,225083+09 

83414. 

49 

3067. 

1.0000 

-1,4180 

.225078+09 

78148, 

bO 

31bO. 

1.0000 

-1,7657 

.225073+09 . 

73280, 

bl 

3213. 

1,0000 

-2.0871 

,225069+09 

68780. 

b2 

327o. 

1.0000 

-2.3841 

.225065+09 

64622. 

53 

3339, 

1.0000 

-2,6587 

.225061+09 

60778. 

b4 

3402. 

1.0000 

-2,9124 

•225057+09 

57226. 

bb 

346b. 

1.0000 

-3.1470 

.225054+09 

53942, 

b6 

3526. 

1,0000 

-3,3639 

.225051+09 . 

50906 « 

b7 

3591. 

1.0000 

-3.5641 

,225048+09 

46102. 

b6 

3654, 

1.0000 

-3,7493 

.225046+09 

45510. 

b9 

3717. 

1,0000 

-3.9206 

.225043+09 

43112. 

t>0 

376o. 

1.0000 

-4,0787 

•225041+09 

40898. 

61 

3843, 

1,0000 

-4.2250 

,225039+09 

38850. 

62 

3906. 

1.0000 

-4.3601 

.225037+09 

36958. 

63 

39b9, 

' 1.0000 

-4.4850 

•225035+09 

35210. 

64 

4032. 

1,0000 

-4.6004 

,225034+09 

33594. 

65 

409o, 

1.0000 

-4,7073 

.225032+09 

32098, 

66 

41b6, 

1.0000 

-4.8059 

.225031+09 

30718. 

o7 

4221, 

i.uooo 

-4,8970 

.225029+09 

29442. 

66 

4264, 

1.0000 

-4,9813 

.225028+09 

28262. 

69 

4347. 

1.0000 

-5,0591 

.225027+09 

27172. 

70 

4410* 

1.0000 

-5.1311 

.225026+09 

26164. 

71 

4473, 

1.0000 

-5.1977 

.225025+09 

25232, 

72 

453o« 

1.0000 

-5,2591 

.225024+09 

24372, 

73 

4599, 

1.0000 

-5,3160 

.225024+09 

23576. 

74 

4662, 

1.0000 

-5, 3666 

.225023+09 

22840. 

75 

472b, 

1.0000 

-5.4171 

,225022+09 

22160. 

76 

478o. 

1,0000 

-5.4620 

.225022+09 

21532, 

77 

48bi. 

1.0000 

-5.5036 

.225021+09 

20950. 

78 

4914, 

1,0000 

-5.5419 

,225020+09 

20414. 

79 

4977. 

1,0000 

-5,5773 

.225020+09 

19918. 

60 

5040. 

1.0000 

-5.6100 

•225019+09 

19460. 
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FhtQUENCY CONTROL SIMULATION 


LOOP GPIN 100. 


STtP 

T 

VM 

V3 

FO 

uelf 


MCSiC 

VOLTS 

VOLTS 

HZ 

. HZ . 

bl 

5103, 

1,0000 

-5,6403 

,225019+09 

19036. 

b 2 

51oo, 

1.0000 

-5.6683 

.22 5019+09 

10644. 

b3 

5 22 k >. 

1,0000 

-5.6941 

.225018+09 

18292, 

b4 

5292, 

1.0000 

-5.7181 

.225018+09 

17946. 

bb 

5355, 

1.0000 

-5,7401 

.225018+09 

17638. 

b6 

5413, 

1,0000 

-5,7606 

.225017+09 

17352, 

b7 

5431, 

1,0000 

-5,7794 

.225017+09 

17088. 

tiQ 

5544 , 

1,0000 

-5,7969 

.225017+09 

16844, 

b9 

5607, 

1,0000 

-5.8130 

•225017+09 

16618. 

90 

567.), 

1,0000 

-5,8280 

.225016+09 

16408, 

91 

5735, 

1.0000 

-5,8417 

.225016+09 

16216. 

92 

579o, 

1.0000 

-5.8544 

.225016+09 

16038. 

9b 

5859, 

1,0000 

-5,8663 

.225016+09 

15872* 

94 

5922. 

1.0000 

-b. 8771 

.225016+09 

15720. 

95 

593b. 

1,0000 

-5,8871 

.225016+09 

15580. 

96 

6046, 

1.0000 

-5,8964 

.22 5015+09 ... 

„„ 15450. 

97 

6111, 

1.0000 

-5.9050 

.225015+09 

15330. 

93 

6174, 

1.0000 

-5,9130 

.225015+09 

15218, 

99 

6237, 

1,0000 

-5.9203 

.225015+09 

15116. 

luO 

6300 • 

1,0000 

-5,9271 

,225015+09 

15020, 

101 

63o3. 

1,0000 

-5,9333 

.225015+09 

14934, 

lu2 

642o , 

1 • 0000 

-5,9391 

.22*015+09 

14852, 

iu3 

3439 , 

1,0000 

-5,944b 

.225015+09 

14/76, 

l04 

6552, 

1,0000 

-5,9494 

.225015+09 

14708. 

10b 

661s, 

1,0000 

-5.954U 

.225015+09 

14644, 

1 U6 

6673, 

1,0000 

-5.9583 

.225015+09 

14584, 

107 

6741, 

1.0000 

-5.9621 

.225015+09 

14530. 

103 

6804 , 

1.0000 

-5,9657 

•225014+09 

14480, 

109 

68o7, 

1.0000 

-5,9691 

.225014+09 

14432. 

110 

b93u * 

1,0000 

-5,9721 

•225014+09 

14390. 

HI 

6993, 

1.0000 

-5.9750 

•225014+09 

14350. 

112 

705o, 

1.0000 

-5,9777 

.225014+09 

14312, 

113 

7119, 

1.0000 

-5.9801 

.225014+09 

14278, 

114 

7162. 

1.0000 

-5.9824 

.225014+09 

14246. 

11b 

7243, 

1,0000 

-5,9844 

,225014+09 

14218. 

116 

7303, 

1,0000 

-5,9964 

.225014+09 

14190, 

117 

73/1. 

1.0000 

-5.9881 

.225014+09 

14166. 

113 

7434, 

1,0000 

-5,9899 

•225014+09 

14142, 

119 

7497. 

1,0000 

-5.9913 

.225014+09 

14122. 

1*0 

75oU • 

1,0000 

-5,9927 

.225014+09 

14102, 



FKt_QUENCY CONTROL SIMULATION 


LOOP GAIN 100, 


STLP 

T 

VM 

V3 

FO 

oelf 


MC5LC 

VOLTS 

VOLTS 

HZ 

HZ 

121 

7623, 

1.0000 

-5.9940 

.225014409 

14084. 

122 

76bo, 

1.0000 

-5,9953 

.225014409 

14066. 

123 

7749, 

1.0000 

-5.9963 

.225014409 

14052, 

124 

7812. 

1.0000 

-5.9974 

.225014409 

14036, 

125 

7875, 

l.UOOO 

-5.9983 

.225014409 

14024* 

126 

7938, 

1.0000 

-5.9991 

.225014409 

14012, 

127 

6001, 

1.0000 

-6.0000 

.225014409 

14000, 

126 

60b4, 

1.0000 

-6.0007 

.225014409 

13990. 

129 

8127. 

1.0000 

-6.0014 

.225014409 

13980. 

130 

8190. 

1.0000 

-6.0020 

.225014409 

13972, 

131 

6253, 

1,0000 

-6,0027 

,225014409 

13962, 

132 

831u, 

1,0000 

-6.0031 

.225014409 

13956, 

133 

8379. 

1.0000 

-6.0037 

.225014409 

13948, 

134 

8442. 

1.0000 

-6,0041 

•225014409 

13942, 

135 

8503. 

1.0000 

-6,0046 

.225014409 

13936, 

136 

85od. 

1.0000 

-6.0050 

.225014409 

13930, 

137 

8631. 

l.UOOO 

-6.0054 

.225014409 

13924, 

13Q 

8694, 

1.0000 

-6,0057 

.225014409 

13920, 

139 

8757. 

1.0000 

-6,0060 

.225014409 

13916, 

140 

862ij. 

1.0000 

-6.0063 

.225014409. 

13912, 

141 

6663. 

1.0000 

-6,0066 

.225014409 

13908, 

142 

694o. 

1.0000 

-6.0069 

•225014409 

13904. 

143 

9009 , 

1,0000 

-6.0070 

.225014409 

13902, 

144 

9072, 

l.UOOO 

-6.0073 

•225014409 

13898. 

14b 

9135. 

l.UOOO 

-6.0074 

.225014409 

13896, 

146 

9198. 

1.0000 

-6.0077 

.225014409 

13892. 

147 

9261. 

1.0000 

-6.0079 

.225014409 

13890. 

148 

9324, 

l.UOOO 

-6.0080 

.225014409 

13888, 

149 

9367. 

1.0000 

-6.0081 

.225014409 

13886. 

150 

9450, 

l.UOOO 

-6,0063 

.225014409 

13884. 

1d1 

9513, 

1.0000 

-6,0084 

.225014409 

13882. 

152 

9575, 

l.UOOO 

-6,0084 

.225014409 

13882, 

153 

9639. 

1.0000 

-6,0086 

.225014409 

13880. 

154 

9702, 

l.UOOO 

-6.0087 

.225014409 

13878. 

155 

9765. 

1.0000 

-6,0087 

.225014409 

13078. 

156 

9826. 

1.0000 

-6.0089 

.225014409 

13876. 

157 

9891. 

1.0000 

-6,0090 

.225014409 

13874. 

158 

9954, 

1.0000 

-6.0090 

.225014409 

13874, 

159 

10017, 

l.UOOO 

-6,0090 

♦225014409 

13874. 

loO 

1U060. 

1.0000 

-6,0091 

.225014409 

13872, 



FREQUENCY control SIMULATION 


LOOP GAIN 100. 


stlp 

T 

VM 

V3 

fo 

UELF 


MC6EC 

VOLTS 

VOLTS 

HZ 

M2 

lol 

10143. 

1,0000 

-6.0091 

.225014+09 

13872, 

lb2 

10206. 

1 ,0000 

-6.0093 

.225014+09 

13870, 

lo3 

10269, 

1.0000 

-6,0093 

.225014+09 

13870. 

1&4 

10332. 

1.0000 

-6,0093 

.225014+09 

13870. 

16b 

10393, 

1.0000 

-6.0094 

.225014+09 

13868. 

166 

10466, 

1,0000 

-6.0094 

.225014+09 

13668. 

lo7 

10521. 

1.0000 

-6.0094 

.225014+09 

13868, 

lod 

10564, 

1.0000 

-6.0094 

.225014+09 

13868, 

169 

10647. 

1.U000 

-6.0096 

.225014+09 

13866. 

170 

10710, 

1,0000 

-6.0096 

,225014+09 

13866, 

171 

10773. 

1,0000 

-6,0096 

,225014+09 

13866. 

172 

lOBoo, 

1.0000 

-6.0096 

.225014+09 

13666. 

1/3 

10899, 

1.0000 

-6.0096 

.225014+09 

13866, 

174 

109u2, 

1,0000 

-6.0096 

.225014+09 

13866. 

175 

11025. 

1.0000 

-6.0097 

.225014+09 

13864, 

176 

11066, 

1,0000 

-6,0097 

.225014+09 

13864, 

177 

11151, 

1,0000 

-6.0097 

.225014+09 

13864, 

178 

11214, 

1.0000 

-6.0097 

.22 5014+09 

13864, 

179 

112/7, 

1.0000 

-6.0097 

.225014+09 

13864. 

160 

11340, 

1,0000 

-6,0097 

.225014+09 

13864, 

161 

11403. 

1.0000 

-6.0097 

.225014+09 

13864. 

162 

11406, 

l.UOOO 

-6.0097 

.225014+09 

13864, 

163 

11529, 

1.0000 

-6.0097 

.225014+09 

13864. 

164 

11592. 

1.0000 

-6.0097 

.225014+09 

13864. 

165 

11655, 

l.uooo 

-6.0097 

.225014+09 

13864. 

166 

11716, 

1.0000 

-6.0097 

,225014+09 

13864. 

167 

11761. 

1.0000 

-6.0099 

.225014+09 

13862, 

166 

11844. 

1.0000 

-6.0099 

.225014+09 

13862. 

169 

11907. 

l.uooo 

°6 • 0099 

.225014+09 

13862. 

190 

11970, 

1,0000 

-6.0099 

.225014+09 

13862. 

191 

12033. 

1.0000 

-6 o 0099 

.225014+09 

13862. 

192 

l209o. 

1,0000 

-6.0099 

•225014+09 

13862. 

193 

12159. 

1.0000 

-6.0099 

•226014+09 

13862. 

194 

12222. 

1.0000 

-6.0099 

• 226014 + 09 

13862. 

195 

12266, 

1,0000 

-6, 0099 

.225014+09 

13662. 

196 

12346, 

1.0000 

-6.0099 

.225014+09 

13062, 

197 

12411, 

1.0000 

-6.0099 

.225°1 4 *09 

13862, 

196 

12474, 

1.0000 

-6.0099 

.225014+09 

13862. 

199 

1253/. 

1.0000 

-6.0099 

•225014+09 

13862, 

2u0 

1260 j , 

1.0000 

-6.0099 

.225014+09 

13862. 


W F Iim 
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APPENDIX F 


FORTRAN program for nonlinear AFC system. 
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STEP 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 
27 
26 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


FREOUENCY CONTROL SIMULATION 



LOOP GAIN 250 * 



T 

VM 

V 3 

FO 

DELF 

MCSEC 

VOLTS 

volts 

HZ 

HZ 

63 . 

• 0000 

- 7.0004 

. 225000+09 

- 2 .. 

126 . 

• 0000 

- 7.0004 

. 225000+09 

- 2 . 

189 . 

.0000 

- 7.0000 

• 225000 + 09 ... 

0 *_ 

252 . 

• 0000 

- 7.0000 

. 225000+09 

0 . 

315 . 

• 0000 

- 7.0000 

. 225000 + 09 -. 

J>*. 

378 . 

• 0000 

- 7.0000 

• 225000+09 . 

Oil. 

441 . 

• 0000 

- 7.0000 

• 225000+09 

0 . 

504 . 

• 0000 

- 7.0000 

. 225000 + 09 . 

0 . 

567 . 

• 0000 

- 7.0000 

. 225000+09 

0 . 

630 . 

1.0000 

12.0000 

• 225800+09 

799998 . 

693 . 

1.0000 

12.0000 

. 225797+09 

797026 .-, 

756 . 

1.0000 

12.0000 

. 225794+09 

793814 . 

819 . 

1.0000 

12.0000 

. 225790 + 09 . 

790362 . 

682 . 

1.0000 

12.0000 

. 225787+09 

786674 . 

945 . 

1.0000 

12.0000 

. 225783+09 

782740 .-, 

1008 . 

1.0000 

12.0000 

. 225779+09 

778572 . 

1071 . 

1.0000 

12.0000 

. 225774+09 

774162 . 

1134 . 

1.0000 

12.0000 

. 225770+09 . 

769514 .- 

1197 . 

1.0000 

12.0000 

. 225765+09 

764632 . 

1260 . 

1.0000 

12.0000 

• 225760 + 09 . 

759512 .-. 

1323 . 

1.0000 

12.0000 

. 225754 + 09 . 

754154 .- 

1386 . 

1.0000 

12.0000 

. 225749+09 

748564 . 

1449 . 

1.0000 

12.0000 

. 225743+09 

742736 .- 

1512 . 

1.0000 

12.0000 

. 225737+09 

736678 . 

1575 . 

1.0000 

12.0000 

© 225730+09 

730380 . 

1638 . 

1.0000 

12.0000 

. 225724+09 

723854 .- 

1701 . 

1.0000 

12.0000 

. 225717+09 

717090 . 

1764 . 

1.0000 

12.0000 

. 225710+09 

710098 *- 

1827 . 

1.0000 

12.0000 

. 225703+09 

702872 . 

1890 . 

1.0000 

12.0000 

. 225695+09 

695414 .- 

1953 . 

1.0000 

12.0000 

• 225688+09 

687724 . 

2016 . 

1.0000 

12.0000 

• 225680+09 

679804 . 

2079 . 

1.0000 

12.0000 

. 225672+09 

671 656 .. 

2142 . 

1.0000 

12.0000 

. 225663+09 

663280 . 

2205 . 

1.0000 

12.0000 

. 225655+09 

654674 .. 

2268 . 

1.0000 

12.0000 

• 225646+09 

645838 . 

2331 . 

1.0000 

12.0000 

. 225637+09 

636778 . 

2394 . 

1.0000 

12.0000 

. 225627+09 

627490 .. 

2457 . 

1.0000 

12.0000 

• 225618+09 

617976 . 

2520 . 

1.0000 

12.0000 

• 225608+09 

608234 .- 
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FREQUENCY CONTROL SIMULATION 
LOOP GAIN 250 * 


STEP 

T 

VM 

V3 


MCSEC 

VOLTS 

VOLTS 

41 

2583 * 

1.0000 

12.0000 

42 

2646 . 

1.0000 

12.0000 

43 

2709 . 

1.0000 

12.0000 

44 

2772 . 

1.0000 

12.0000 

45 . 

2835 . 

1.0000 

12.0000 

46 

2898 * 

1.0000 

12.0000 

47 

2961 . 

1.0000 

12.0000 

48 

3024 . 

1.0000 

12.0000 

49 

3087 . 

1.0000 

12.0000 

50 

3150 . 

1.0000 

12.0000 

51 

3213 * 

1.0000 

12.0000 

52 

3276 . 

1.0000 

12.0000 

53 

3339 . 

1.0000 

12.0000 

54 

3402 . 

1.0000 

12.0000 

55 

3465 . 

1.0000 

12.0000 

56 

3528 . 

1.0000 

12.0000 

57 

3591 . 

1.0000 

12.0000 

58 

3654 . 

1.0000 

12.0000 

59 

3717 . 

1.0000 

12.0000 

60 

3780 . 

1.0000 

12.0000 

61 

3B43 * 

1.0000 

12.0000 

62 

3906 * 

1.0000 

12.0000 

63 

3969 . 

1.0000 

12.0000 

64 

4032 . 

1.0000 

12.0000 

65 

4095 . 

1.0000 

12.0000 

66 

4158 . 

1.0000 

12.0000 

67 

4221 . 

1.0000 

12.0000 

68 

4284 . 

1.0000 

12.0000 

69 

4347 . 

1.0000 

12.0000 

70 

4410 . 

1.0000 

12.0000 

71 

4473 . 

1.0000 

12.0000 

72 

4536 . 

1.0000 

12.0000 

73 

4599 . 

1.0000 

12.0000 

74 

4662 . 

1.0000 

12.0000 

75 

4725 . 

1.0000 

12.0000 

76 

4788 . 

1.0000 

12.0000 

77 

4851 . 

1.0000 

9.7018 

78 

4914 . 

1.0000 

6.9054 

79 

4977 . 

1.0000 

4.5846 

60 

5040 . 

1.0000 

2.6661 


FO DELF 

HZ HZ 

• 225598+09 ... 598266 *.. 

• 225588+09 588072 * 

• 225578+09 577658 * 

. 225567+09 567016 * 

• 225556+09 556156 * 

• 225595+09 595066 * 

• 225539+09 533756 * 

• 225522+09 522229 * 

• 225510+09 510970 . 

* 225998+09 498996 . 

. 225986+09 986300 * 

• 225979+09 973892 . 

* 225961+09 961256 * 

• 225498+09 998402 * 

. 225935+09 435332 * 

. 225422+09 422042 * 

• 225409+09 408534 * 

. .. . 225395+09 394812 . 

• 225381+09 380874 * 

• 225367 + 09 .... 366718 * 

. 225352+09 352348 * 

. 225338+09 337764 * 

. 225323+09 322962 * 

* 225308+09 307948 * 

. 225293+09 292724 * 

. 225277+09 277282 * 

• 225262+09 261630 * 

. 225246+09 245764 . 

• 225230+09 229690 * 

. 225213+09 213400 . 

. 225197+09 196904 . 

* 225180+09 180198 * 

. 225163+09 163280 * 

• 225146+09 146156 * 

. 225129+09 126822 * 

* 225111+09 111280 * 

. 225094+09 93 $ TOr “~ 

. 225078+09 77870 . 

• 225065+09 64874 . 

. 225054+09 54130 . 
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FREQUENCY CONTROL SIMULATION 


LOOP GAIN 250* 


STEP T VM V3 FO DELF 

MCSEC VOLTS VOLTS HZ — HZ „ 

81 5103* 1.0000 1.0846 . 225045*09 45274*. 

82 5166. 1.0000 -.2150 .225038*09 37996. 

83 5229. 1.0000 -1.2807 .225032*09 32028*_ 

04 5292. 1.0000 -2.1536 .225027*09 27140. 

85 5355. 1.0000 -2.8671 .225023*09 .23144*. 

66 5418. 1.0000 -3.4493 .225020*09 19884*.- 

87 5481. 1.0000 -3.9246 .225017*09 17222* 

68 5544. 1.0000 -4.3118 .225015*09 15054*. 

89 5607. 1.0000 -4.6279 .225013*09 132B4* 

90 5670. 1.0000 -4.8846 *225012*09 .11846*.. 

91 5733. 1.0000 -5.0939 .225011*09 — - 10674*- 

92 5796. 1.0000 -5.2643 .225010*09 9720. 

93 5859. 1.0000 -5.4025 .. .225009*09 6946*- 

94 5922. 1.0000 -5.5150 .225008*09 8316. 

95 5985. 1.0000 -5.6064 .225008*09- „..7804*-- 

96 6048. 1.0000 -5.6811 .225007*09 7386* 

97 6111. 1.0000 -5.7414 .225007*09 7048. 

98 6174. 1.0000 -5.7907 .225007*09 6772*_ 

99 6237. 1.0000 -5.8307 .225007*09 6548* 

100 6300. 1.0000 -5.8629 . .225006*09 6368. 

101 6363. 1.0000 -5. 8896 .225006*09 6218. 

102 6426. 1.0000 -5.9111 .225006*09 6098. 

103 6489. 1.0000 -5.9282 . .225006*09 6002*- 

104 6552. 1.0000 -5.9425 .225006*09 5922* 

105 6615. 1.0000 -5.9543 .225006*09 5856.- 

106 6678. 1.0000 -5.9636 .225006*09- 5804*- 

107 6741. 1.0000 -5.9711 .225006*09 5762. 

108 6804. 1.0000 -5.9775 .225006*09 5726. _ 

109 6867. 1.0000 -5.9821 .225006*09 5700. 

110 6930. 1.0000 -5.9861 .225006*09 5678*. 

111 6993. 1.0000 -5. 9896 .225006*09 5658* 

112 7056. 1.0000 -5.9925 .225006*09 5642* 

113 7119. 1.0000 -5.9943 .225006*09 5632*- 

114 7182. 1.0000 -5.9961 .225006*09 5622. 

115 7245. 1.0000 -5.9979 .225006*09 5612*. 

116 7308. 1.0000 -5.9989 .225006*09 ... - 5606* 

117 7371. 1.0000 -6.0000 .225006*09 5600* 

118 7434. 1.0000 -6.0004 .225006*09 .. 5598. 

119 7497. 1.0000 -6.0014 .225006*09 5592* 

120 7560. 1.0000 -6.0014 .225006*09 5592* 
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APPENDIX G 


FORTRAN program for AFC system with alternate 
sampler with limiter and nonlinear VCO model. 
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ENo of OMplLATiON: no oIAGNOSTicS 



Loop g in 


250 


STfcP 

T 

VM 

V3 

PO 

DFLF 


MCSEC 

VOlTS 

VOLTS 

HZ 

HZ 

0 

-63.49 

• oOOO 

-7,0004 

.225^0+09 

-2* 

0 

-63.34 

.0000 

-7,0004 

• 225000+09 

-2. 

0 

-63.19 

.0000 

-7,0004 

,225000+09 

-2. 

0 

-63.04 

.0000 

-7,0004 

.225000+09 

-2. 

0 

-62.89 

.0000 

-7,0004 

,225000+09 

-2* 

0 

-62.74 

• oooo 

-7,0004 

•225000+09 

-2. 

0 

-62.59 

. OOOO 

-7,0004 

• 2250*0+09 

-2. 

0 

-62.44 

• OOOO 

-7,0004 

•225000+09 

-2, 

0 

-62.29 

.0000 

-7,0004 

•2250n0+09 

-2* 

0 

-62.14 

• oooo 

-7,0004 

.2250n0+09 

-2. 

0 

-61.99 

.0000 

-7,0004 

.225000+09 

-2* 

1 

. oO 

1.0000 

12,0000 

•2258n0+09 

7o99«8* 

1 

• lb 

1 * oooo 

12,0000 

•2258n0+09 

7o97i2, 

1 

• 30 

l.oono 

12.0000 

•225799+n9 

79*422, 

1 

.4b 

1*0000 

12.0000 

•225799+09 

79*132. 

1 

• 60 

1.0000 

12,0000 

•225799+09 

708638# 

1 

.75 

1.0000 

12,0000 

.225799+09 

798544 * 

1 

• 90 

1.0000 

12,0000 

,225798+09 

798248# 

1 

1.05 

1.0000 

12,0000 

•225 7 98+09 

797946# 

1 

1.20 

1.0000 

12,0000 

. 225798+09 

797640# 

1 

1.35 

1.0000 

12,0000 

.225797+09 

797336# 

1 

1.50 

1.0000 

12,0000 

,225797+09 

797026* 

2 

63.49 

1.0000 

12,0000 

.225797+09 

797026* 

2 

63 . 64 

1.0000 

12,0000 

.225797+09 

796718. 

2 

63.79 

1.0000 

12,0000 

.225796+09 

796404* 

2 

63.94 

1.0000 

12,0000 

.225796+09 

796092* 

2 

64.u9 

1.0000 

12,0000 

.225796+09 

79*5772* 

2 

64.24 

1*0000 

12,0000 

.225795+09 

795452. 

2 

64.39 

1.0000 

12,0000 

•225795+09 

795130. 

2 

64.54 

1.0000 

12,0000 

.225795+09 

794802* 

2 

64.69 

1.0000 

12,0000 

.225794+09 

794474* 

2 

64.84 

1.0000 

12,0000 

•225794+09 

794146. 

2 

64.99 

1.0000 

12,0000 

.225^94+09 

793814. 

3 

126.98 

1.0000 

12,0000 

.225794+09 

793814* 

3 

127.13 

1.0000 

12 , OOOO 

, 225" r 93+09 

7q34*2« 

3 

127,28 

1.0000 

12,0000 

•225 7 93+09 

793148. 

3 

127,43 

1.0000 

12,0000 

• 225793+09 

792804. 

3 

127.58 

1.0000 

12.0000 

.225^92+09 

792460* 

3 

127.73 

1.0000 

12,0000 

.225792+09 

792118# 

3 

127. 88 

1.0000 

12,0000 

.225^92+09 

791774* 

3 

1 2® • 0 3 

1.0000 

12,0000 

,225 7 9l+09 

791426# 

3 

128,18 

1.0000 

12.0000 

.225791+09 

791072. 

3 

128,33 

i.0000 

12.0000 

.225791+09 

790718. 

3 

12&.48 

1.0000 

12,0000 

.225790+09 

790366# 
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loop gain 250 


STEP 

T 

VM 

V3 

FO 

delf 


MCSEC 

VOlTS 

VOLTS 

HZ 

HZ 

4 

190*46 

1 

0000 

12.0000 

.225790+09 

7Q0366. 

4 

190*63 

1 

0000 

12,0000 

.225790+09 

790004* 

4 

190*78 

1 

0000 

12,0000 

.225790+09 

709646* 

4 

190.93 

1 

0000 

12.0000 

.225789+09 

7ft9282. 

4 

191.06 

1 

oono 

12,0000 

.225789+09 

7089*6* 

4 

191.23 

1 

0000 

12,0000 

•225789+09 

708550. 

4 

191.38 

1 

0000 

12,0000 

•225788+09 

700176* 

4 

191.53 

1 

0000 

12,0000 

.225788+09 

707802* 

4 

191.66 

1 

0000 

12,0000 

.225787+09 

787420. 

4 

191.63 

1 

0000 

12,0000 

•225787+09 

707052. 

4 

191.96 

1 

0000 

12,0000 

•225787+09 

706670* 

5 

253,97 

1 

0000 

12,0000 

.225787+09 

706670* 

5 

2&4*i2 

1 

0000 

12.0000 

•225766+09 

706292. 

5 

2 5 4.27 

1 

0000 

12,0000 

•225786+09 

705902. 

5 

25^.42 

1 

0000 

12,0000 

.225786+09 

705510. 

5 

254.57 

1 

0000 

12.0000 

•225785+09 

705120. 

5 

254.72 

1 

0000 

12,0000 

•225785+09 

704736. 

5 

254.87 

1 

0000 

12.0000 

.225784+09 

704342. 

5 

255.o2 

1 

0000 

12,0000 

•225784+09 

703946* 

5 

255 , \1 

1 

0000 

12.0000 

•225784+09 

703546* 

5 

255.32 

1 

oooo 

12,0000 

•225783+09 

703144* 

5 

255,47 

1 

0000 

12,0000 

•225783+09 

702740. 

6 

3l7.4b 

1 

0000 

12,0000 

.225783+09 

702740* 

6 

3l7,61 

1 

oooo 

12,0000 

.225782+09 

702336* 

6 

317.76 

1 

0000 

12,0000 

.225782+09 

791924. 

6 

3i7,91 

1 

oooo 

12.0000 

.225782+09 

701512* 

6 

3i8,ob 

1 

ooco 

12,0000 

,225781+09 

701100* 

6 

3l8,2l 

1 

oooo 

12,0000 

•225781+09 

700604. 

6 

318,36 

1 

oOOO 

12,0000 

.225780+09 

700264. 

6 

3i8,51 

1 

oooo 

12,0000 

•225780+09 

77901*6. 

6 

3i8.66 

1 

oooo 

12,0000 

#225779+09 

779422* 

b 

318.Q1 

1 

oooo 

12,0000 

•225779+09 

778990. 

b 

3i8.96 

1 

oooo 

12,0000 

•225779+09 

770572* 

7 

380,95 

1 

oooo 

12,0000 

•225779+09 

770572. 

7 

381.10 

1 

OOPO 

12.0000 

.225778+09 

770140* 

7 

381,25 

1 

oooo 

12,0000 

.225778+09 

7777n6. 

7 

381.40 

1 

oono 

12,0000 

.225777+09 

777270. 

7 

381.55 

1 

oooo 

12.0000 

.225777+09 

776036* 

7 

3fil * 70 

1 

oooo 

12,0000 

.225776+09 

776396. 

7 

381.85 

1 

OOOO 

12,0000 

.225776+09 

775954. 

7 

3e2 . jO 

1 

0000 

12,0000 

•225776+09 

775500. 

7 

302,15 

1 

oOOO 

12,0000 

•225775+09 

775066. 

7 

382.30 

1 

oooo 

12,0000 

•225775+09 

774612. 

7 

382.45 

1 

oooo 

12,0000 

.225774+09 

774150. 
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loop .a in 250 


STtP 

T 


VM 

V3 

FO 

oflf 


MCSEC 

VOLTS 

VOLTS 

HZ 

HZ 

a 

445.94 

1 

• 0000 

12.0000 

.225770+09 

7**514. 

9 

5()9.44 

1 

.0000 

12,0000 

.225765+09 

7*4628. 

10 

572.93 

1 

.0000 

12.0000 

•225760+09 

75*510. 

11 

636,42 

1 

• 0000 

12,0000 

,225754+09 

75*158. 

12 

699.91 

1 

• 0000 

12.0000 

,225749+09 

74*560. 

13 

763,40 

1 

• oono 

12.0000 

.225743+09 

742756. 

14 

826,90 

1 

• 0000 

12,0000 

.225737+09 

736674. 

15 

890.39 

1 

• 0000 

12,0000 

.225730+09 

730376. 

16 

953,08 

1 

• 0000 

12,0000 

♦225724+09 

723854. 

17 

1017,37 

1 

• 0000 

12.0000 

♦225717+09 

717090. 

18 

1080.07 

1 

• 0000 

12,0000 

.225710+09 

710094. 

19 

1144.36 

1 

♦ oooo 

12.0000 

.2257-3+09 

7*28*8. 


12o7.05 

1 

• OOOO 

12.0000 

.225695+09 

69541ft, 

21 

1271.34 

1 

.0000 

12,0000 

.225688+09 

6077?4. 

22 

1334.03 

1 

.0000 

12,0000 

♦225680+09 

679804. 

23 

1398.33 

1 

• OOOO 

12,0000 

.225672+09 

671656. 

24 

1461.82 

1 

.0000 

12,0000 

.225663+09 

6*3280* 

25 

1525,31 

1 

.0000 

12.0000 

.225655+09 

654674. 

26 

1508,00 

1 

.0000 

12.0000 

.225646+09 

645838. 

27 

1652,29 

1 

• oooo 

12,0000 

.225637+09 

636774. 

2a 

17i5.79 

1 

• nooo 

12.0000 

.22562^+09 

627406. 

29 

1779.28 

1 

.0000 

12.0000 

.225619+09 

617972. 

3q 

1842.77 

1 

.0000 

12.0000 

.2256C8+09 

6ft8230. 

31 

19o6.26 

1 

• oooo 

12.0000 

•225598+09 

598262. 

32 

l9 6 9./5 

1 

.0000 

12.0000 

.225588+09 

548068. 

33 

2033,25 

1 

• oooo 

12.0000 

•225578+09 

577650. 

34 

2096,74 

1 

.0000 

12.0000 

.225567+09 

5*7oi2. 

35 

21bU.23 

1 

,0000 

12,0000 

.225556+09 

556148* 

36 

2223,72 

1 

• oooo 

12,0000 

.225545+09 

545058. 

37 

2207,21 

1 

.0000 

12.0000 

.225534+09 

533752. 

38 

2350.71 

1 

• oooo 

12.0000 

.225522+09 

522216. 

39 

2414,20 

1 

• oooo 

12.0000 

•225510+09 

510462, 

40 

2477.69 

1 

.0000 

12,0000 

.225498+09 

408492, 

41 

2541.18 

1 

.0000 

12.0000 

.225486+09 

406296, 

42 

26q4,67 

1 

• oooo 

12,0000 

.225474+09 

473880. 

43 

2668,17 

1 

• oooo 

12,0000 

•225461+09 

461250. 

44 

27j1,66 

1 

.0000 

12,0000 

.225448+09 

448394. 

45 

2795.15 

1 

• oooo 

12.0000 

•225435+09 

435320. 

46 

2858.64 

1 

.0000 

12.0000 

•225422+09 

422034. 

47 

29 2 2.13 

1 

.0000 

12,0000 

•2254Q9+09 

40*526. 
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Loop gain 250 


STfcP 

T 

VM 

V3 

FO 

OIL ? 


MCSEC 

VOlTS 

VOLTS 

HZ 

HZ 

48 

29e5.63 

1.0000 

12.0000 

.225395+09 

394804* 

49 

3049.12 

1.0000 

12,0000 

•225381+09 

3n0862* 

50 

3112.61 

1.0000 

12,0000 

.225367+09 

366710* 

51 

3176,10 

1.0000 

12,0000 

.225352+09 

352336* 

52 

3239.60 

1.0000 

12.0000 

•225338+09 

337752. 

53 

33o3. 09 

1.0000 

12,0000 

.225323+09 

3*2952* 

54 

3366.56 

1.0000 

12,0000 

.225308+09 

3^7936. 

55 

3430.07 

1.0000 

12,0000 

.225293+09 

292708* 

56 

3493.56 

1.0000 

12.0000 

•225277+09 

277270* 

57 

3557 . 06 

1.0000 

12,0000 

•225262+09 

261614* 

58 

3620.55 

1.0000 

12,0000 

,225246+09 

245750. 

59 

3684.04 

1.0000 

12.0000 

.225230+09 

229678. 

50 

3747.53 

1.0000 

12.0000 

•225213+09 

213386* 

tol 

3011. (j2 

1.0000 

12.0000 

•225197+09 

196888* 

52 

3874.52 

1.0000 

12.0000 

•225180+09 

180182* 

53 

3938,ol 

1.0000 

12.0000 

•225163+09 

163264* 

54 

4001.50 

1.0000 

12,0000 

•225146+09 

146140* 

55 

4064.99 

1.0000 

12.0000 

.225129+09 

120804* 

56 

4128.46 

1.0000 

12.0000 

• 225U1 + H9 

111262. 

57 

4191.96 

1.0000 

9,8221 

•225094+09 

94204. 

58 

4255.47 

1.0000 

7.2204 

.225080+09 

79634* 

69 

4316,96 

1.0000 

5.0282 

.225067+09 

67358. 

70 

4302,45 

1.0000 

3,1868 

•225057+09 

57046. 

71 

4445,94 

1.0000 

1,6439 

.225048+09 

4*406* 

72 

45o9,44 

1.000O 

,3539 

.225041+09 

41102* 

73 

4572.93 

1.0000 

-.7225 

•225035+09 

30154* 

74 

4636,42 

1.0000 

-1,6196 

.225030+09 

30130. 

75 

4699.91 

1.0000 

-2,3666 

•225026+09 

25946* 

76 

47e3,40 

1.0000 

-2.9871 

•225022+09 

22472. 

77 

4826,90 

1.0000 

-3,5032 

•225020+09 

19582. 

78 

4890.39 

1.0000 

-3.9314 

.225017+09 

17184. 

79 

4963.86 

1.0000 

-4,2871 

.225015+09 

15192. 

80 

5017.37 

l.OOpO 

-4.5814 

•225014+09 

13544* 

51 

5080.86 

1.0000 

-4,8257 

•225012+09 

12176* 

52 

5144,36 

1.0000 

-5.0286 

•225011+09 

11040. 

53 

52o7.85 

1.0000 

-5,1964 

.225010+09 

lnioo. 

64 

5271.34 

1.0000 

-5,3354 

•22500^+09 

9322. 

55 

5334.93 

1.0000 

-5.4500 

•225009+09 

0680* 

66 

5398.33 

1.0000 

-5,5454 

•225008+09 

0146* 

57 

5461.32 

1.0000 

-5,6246 

•225008+09 

7702. 
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loop 6 in 250 , 


STtP 

T 

VM 

V3 

PO 

OELP 


mcsec 

VOlTS 

VOLTS 

HZ 

HZ 

88 

55^5,31 

1*0000 

-5,69n 0 

•2250p7+09 

7356. 

89 

556ti.au 

1 *0000 

-5,7439 

.225007+09 

7os*. 

*0 

5652*29 

1.0000 

-5.7866 

•225007+09 


*i 

5715.79 

1.0000 

-5.8254 

• 2250f?7+09 

6578, 

92 

5779.28 

1.0000 

-5,8564 

.225006+09 

6*e*. 

93 

5842.77 

1.0000 

-5,8818 

•225006+09 

6262, 

94 

59n6.26 

1.0000 

-5,9029 

* 2250n6+09 

6!**, 

95 

5969.75 

1*0000 

-5,9200 

•225006+09 

60 * 8 . 

96 

6033.25 

l.oOCQ 

-5.934 6 

•225006+09 

5966. 

97 

6096.74 

1.0000 

-5.9464 

•225006+09 

5900. 

98 

6160*23 

l.ooro 

-5,9561 

*2250n6+09 

58*6. 

99 

6223.72 

1.0000 

-5. 964 6 

•225006+09 

579., 

luo 

6287.21 

1.0000 

-5,9711 

•225006+09 

5762, 
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APPENDIX H 

LISTINGS OF ALSAP ROUTINES 

In addition to the ALSAP main routine, there are special 
subroutines for numerical integration (NINTEG) , transfor- 
mation matrix generation (GREG) , and loading the input 
data (DATA IN ) . These routines are listed in this appendix. 
Also used but not listed here are two utility routines for 
finding the eigenvalues and eigenvectors of a matrix (MEGVEL 
and IGVEC) . 
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COMMENT CCS BANDPASS FILTER S/12/71 SKk kEShunSE 

comment .*#«#»« may 25*1971 ••*«*«•*»• 

C APPROXIMATE linear jXSlEM ANALYSIS PROGRAM (AL^AP) 

COMMON A(13.13).o(13,13). NuR I G *M * N * a \ 1 3 ) • NE v i PR • C i 1 3 t 1 3 ) * 0 * 1 3 * 1 J > • 

1 LAMWEL*LAMIMG*SCALE fPVI 13) fNLt'asl #1 lEMP.^i 131 .IPOPi . IPR.FG 
R^AL L AMPFL » LAM | MG 
DOUBLE PRECISION PV 

C*************" . 

DIMENSION RSR(13t *R3I ( 13) *CMAT( 1 3*13 ) *KSKNU3)*KSlN(i 3 I 

DOUcLE PRECISION VR( 13) *Vl ( 13 ) *RR*RI • lEMPD 

CALL DATAIN 

NOR I G * N 

DO S011 I « 1 * N 

DO SO 10 J a l.N 

5010 C( I * J ) * 0. 

5011 C ( 1 • I ) « 1* 

DO 5012 I * 1 * N 
DO 5012 J * 1 *M 

5012 D ( I ♦ J > = 0. 

C * * * FIND THE LARGEST AND SMALLEST NON-^ERO MAGNIluOE EmEMENis OF MAiRIx A 
C FOR DETERMINING SCALING 10 BE USED 10 JMPROvE ‘i-jMEPlCAL ACCuRACi 

C IN THE simplification PROCESS. 

TEMPI = 0. 

TEMP 2 * l.E+30 
DO 6000 I = 1 • N 
DO 6000 J = 1 * N 
TEMP = ABS ( A( I * J ) ) 

IF (TEMP.LT . 1 .E-25 ) GO TO 6000 
IF( TEMP. GT. TEMPI > TEMPI a TEMP 

IF (TEMP.LT* T£MP2 ) TEMP2 = iEmp 

6000 CONTINUE 

I = AL0G10< TEMPI ) 

J = ALOGlOt TEMP2 ) 

K = I + ( I+J) /2. 

SCALE = ( 1 0*0 ) **K 
WRITE (6.6) 

WR I TE ( 6 .6001 ) SCALE 

6001 FORMAT (10X.48HTHE PROGRAM HAS CHOSEN THE SCALE FACIOR 10 BE * *E12 
1 .4 ) 

CO M MENT*******PR I NT WARNING RELATIVE 10 SCALE FAC»u* AnD E T GEn v AluEo**** 

WR I TE I 6*6002 ) 

WR I TE I 6 *6003 1 
WR I TE ( 6*6004 > 

WRI TE ( 6.6005 ) 

6002 FORMAT ( 10X.51H*#*N0TE— ALL PRINTED VALUES Fuk EIGENVALUES AnD F{JK ) 

6C03 FORMAT ( 1CX *51H***ELEMENTS OF THE nATnlX A n,uST 3E mUlUklIED By) 

6004 FORMAT ( 10X *51H*#*THE ABOVE SCAlE F AC Tun BEFunE E X TEn^Ai. USE. »-*-») 

6CG5 FORMAT ( 10X *54h**»N0 OTHER OLANT I T I ES IN THE PkInTjU I kEoUIkE SCAlI 

1NG) 

COMMENT ***** END OF PRINT WARNING****** 

WRI TE ( 6.6 ) 

WRITE ( 6*7000) 

7000 FORMAT ( 1H1 ) 

SCAL1V = 1. /SCALE 
DO 40 I a i ,n 
DO 40 J a 1 ,N 
40 A ( I , J) = All . J)*SCALIV 
WRI TE (6.6 ) 

DO 999 IROOT a 1 * NEVT BR 
IFUPRTFG.EO.O ) GO TO 8001 
WR I TE ( 6 .6 ) 

WRITE* 6.6006 ) 

6006 FORMAT ( 10x .31HEIGENVALUES OF PRESENT MATRIX A) 

8001 CONTINUE 
WR I TE < 6 *6 1 

CALL MEGVAL* N.A.RSR.RSI ) 

WRITE (6.6) 

TEMP a RSR(l) 

I TEMP * 1 
DO 1100 I = 2 * N 
I F ( TEMP • LT • RSR(I)) GO TO 1100 
TEMP = R SR ( I ) 

ITEMP = I 
1100 CONTINUE 

LAMREL = RSR( I TEMP )*SCAlE 
LAMIMG a RSI ( I TEMP )*SCALE 
RR a R SR ( I TEMP ) 

R I a RSI ( ITEMP ) 

CALL I GVEC ( N * 1 3 .A .RR ,R I ,VK *VI ) 

I FI I PRTFG.EQ.O) GO TO 8002 
WRITEI6.42 ) 

42 FORMAT I 10X . 38HE IGENVALUE WITH SMALLEST TImE CunSIAnI) 

WRITEI6.2) ITEMP, KR.RI 

WRI TE ( 6 .6 ) 

WRI TE < 6.43 ) 

43 FORMAT ( 10X .41HE IGtNVECTOK FUm THE ABuVE FAST E!GE*VA L UE) 

WRI TE (6 .6 ) 

DOS J = l.N 
5 WRITE (6.2) VRUliVlIJ) 

8002 CONTINUF 
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CG , wr\T--iva i tc a * v for cc'car Ison with lambda ■ v 
h.gvfc = 4 

6500 CONTINUE 

•a H ! TE ! 6 *0 I 

WPI Tf( 1 E OvE C * 44 ) 

44 F OPM AT ( 1 OX .25hE I GENVEC TOR VALIDITY TEST) 

WRJ TE ( 6 .6 ) 

DO 32 = 1 • N 

PRJNTk = pR«VR(pC*) - H*V1KM 
PRINT1 = RR*VI(K<) ♦ RI#VR(<IC) 

VLAMR = 0. 

DO 30 < = l.N 

30 VLAV.R = VLAMR * AMC<«K )«VR(<) 

vlami = o, 

DO 31 < = l.N 

31 .VLAMI = VLAMI + A ( K< * < ) *V I CM 

I F ( 1 ABS! VLAVR-PRINTR > .GT.ABS! VlAmk*1.E-3) WAmD.I IEGVFC.nE. 6) ) Gu 
1 TO 6501 

I FI (Abb I VLAMI -PRINT I ) *GT.AB$ ( VLAMI *1 *£-3 >) .AM). ( I EGVFC .nE.6 ) ) GU 
1 TO 6501 
GO TO 6502 

6501 IEGVEC = 6 
GO TO 6500 

6502 CONTINUE 

WRITE! IEGVEC. 2 ) VLAMR ,PR INTR.K<, VLAMI .PRINT I 

32 CONTINUE 
WRITE (6. 6) 

TEMP = DABS ( VR ( 1 ) ) 

I TEMP = 1 
DO 1101 K = 1 • N 

IF! TEMP.GT.DA6$(VR{<) ) ) GO TO 1101 
TEMP = DABSJVR (<) ) 

I TEMP = K 

1101 CONTINUE 
TEMPD = VR ( N ) 

VR(N) = VR! ITEMP) 

VR(ITEMP) = TEMPD 

IF ( IPRTFG.EO.O ) GO TO 0003 
wR I T E ( 6 .49 ) ITEMP. N 

49 FORMAT { 10X .20HROWS (AND COLUMNS) .12. 7H A<nD .I2.36H xjf wAIkIx 
1A HAVE BEEN INTERCHANGED ) 

8003 CONTINUE 

COMMENT INTERCHANGE ROWS AND COLS. OF MATRIX A 

DO 1102 K l.N 
TEMP * A(K.N) 

A ( < * N ) = A ( K * I TEMP ) 

1102 AU. ITEMP) = TEMP 
DO 1103 K - l.N 
TEMP = A(N*K) 

A ! N * <. ) = A ( ITEMP, <) 

1103 A ( I TEMP * K ) = TEMP 
DO 1200 J = l.M 
TEMP = B ( N * J ) 

6 1 N * J ) = 3 t I TEMP * J ) 

1200 31 I TEMP, J) = TEMP 
WRITE(6.6) 

NLESS 1 = N - 1 
WRITE (6. 6) 

I F ( IPRTFG.EO.O ) GO TO 0004 
WRITE (6 *45 ) 

45 FORMAT ( 10x .A6HFIRST N-l COMPONENTS OF NORMALIZED E I GFNVECTOR ) 

0004 CONTINUE 

DO 1001 I = 1.NLESS1 

1001 PV ( I ) = VR ( 1) /VR ( N ) 

IF( IPRTFG.EO.O ) GO TO 0005 

w R I T E ( 6 • 2 ) (PV(J).J = 1 .NLESS 1 ) 

WRITE (6.6) 

0005 CONTINUE 

I F I ABS ( LAMIMG ) .L T .A6S ( LAMR£L*1 .E-25 ) ) Gu To 5000 
wRITE (6.5015) 

5C15 FORMAT ( 1 OX .26HC0MPLEX ROOT To 5E KEv.uVED ) 

GO TO 998 
5000 CALL GREG 

DO 1000 I = 1 * NLESS 1 
A! I .N) = 0. 

DO 1000 J = 1.NLESS1 
1000 4(1, J) = A(I.J) - PV!I)*A(N,J) 

A(N.N) = LAMREL/SCALE 
IF! IPRTFG.EO.O ) GO TO 8006 
wRITE 16*46) 

46 FORMAT { 10X .20HTRANSFCRMED MATRIX A) 

WRITE (6.6) 

DO 1002 I - l.N 

1002 w R 1 T F (6*2) (A(I.J).J = l.N) 

w r i r r ( 6 * 6 > 

m R I T L' ( 6 . 4 7 ) 

47 FORMAT ( 10X.35HE IGLNVALUES OF TRANSFORMED MATRIX A) 

CALL MKiVAUN. A *RSR *RS1 ) 

00 J6 CONTINUE 
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COMMENT — GENERATE (N-l) MY (N-l) MATRIX AND FIND ITS EIGENVALUES 
DO 1205 I - ltNLESSl 
DO 1 205 J = 1 * NLLSS1 
CM AT (I*Ji = A ( I . J > 

WR I TE < 6*6 I 

DO 1203 I * 1 t MESS 1 
DO 1203 J = l.M 
b(I *JI = B (I * J) - PV( I ) * 3 < N » J ) 

IF ( IPRTFG.cO.O ) GO TO 8007 
wRI TE (6*6 ) 

WRITE (6*1204) 

FORMAT { 10X.29HMATRIX 6 AFTER T R ANSFORMAT I ON ) 
wRI TE (6.6 ) 

DO 1207 I = 1#N 
WRITE<6.2> ( B( I • J ) * J * 1 * M ) 

WR I T E ( 6 * 1500 ) 

FORMAT ( 1 H 1 J 
CONTINUE 
N - N - 1 
DO 1502 I * 1 ♦ N 
DO 1 502 J = l.N 
AlIiJ) = CMAT ( I * J > 

CONTINUE 
WRITE (6.6) 

*RI TE (6*7000) 

WRI TE ( 6*1206) 

FORMAT < 10X *29HEIGENVALUES CF FINAL MATRIX A) 

CALL MEGVAL ( NL ESS 1 * CMAT *R$RN»RS I N ) 

COM.VENT******THE MATRICES NEEDED FOR DIFFERENTIAL EQUATION SOLUTION 


C********#******#*****aRE PREPARED HERE************ 

C THE MATRIX A MUST 6E UN-SCALED**** 


DO 3000 I = 1*NLESS1 
CO 3000 J = 1 * NLESS 1 
3000 A I I * J ) s A ( I . J ) *SCALE 
C******THE MATRIX 6 IS CORRECT-AS IS******* 

C ******* thE MATRICES C AND C ARE USED TO RECOVER THE ORIGINAL STATE FROM 

C#**** THE SOLUTION OF THE REDUCED DIFFERENTIAL EOUATIUN 

C* * * * (X = C *Z +D*U ) * ****** 

COMMENT* * * INITIALIZATION ****** 

DO 3004 I - 1 »N 

3004 Z ( I ) = 0. 

C* * * * PRINT MATRICES* ****** # 

WRI TE ( 6*6 ) 

WRITE (6.3005) 

3005 FORMAT ( 10X .49HMATR ICES FOR NUMERICAL INTEGRATION AND CONVERSION) 
WRI TE ( 6 . 3006 ) 

3006 FORMAT ( 10X . 19HZD0T=AZ+BU. X=CZ+DU) 

WRI TE < 6 .6 ) 

WRITE (6.6) 

WRITE (6.3007) 

3007 FORMAT ( 10X .8HMATRIX A) 

DO 3008 I = l.N 

3008 WRITEI6.2) (AU.J).J = 1*N) 

WR I T E ( 6 . 6 ) 

WRI TE<6. 3009) 

3009 FORMAT I 1 OX .8HMATRI X B) 

DO 3010 I = l.N 

3010 W R I T E ( 6 * 2 ) ( 3 ( I . J ) * J= 1 *M ) 

WR I TE ( 6 *6 ) 

WRI TE(6.301 1 ) 

3011 FORMAT ( 10X .8HMATRIX C) 

DO 3012 I = 1. NORIG 

3012 W R I T E ( 6 * 2 ) (C(I.J).J = l.N) 

WRI TE ( 6.6) 

WRI TE (6*3013 ) 

3 C 1 3 FORMAT ( 10X .8HMATRIX D) 

DO 3014 I = 1 .NORIG 
3014 WRITEI6.2) ( D( I . J ) * J= l.M) 

C * * * * *END CF MATRIX PRINTING* * * * 

2 FORMAT ( ) 

6 FORMAT (1H ) 

call N I NTEG I N* NOR IG.M.A.B.C .D.Z ) 

998 CONTINUE 
STOP 
FND 

SUBROUTINE N 1 NT EG ( N . NOR IG . V . 4 . B * C *D .2 ) 

DIMENSION A(13, 13). 6(13. 13). 0(13*13) .0(13.13) »Z( 1 3 ) .X ( 1 3 ) .2D0T ( 13 ) 
l.U(l?l 

C DIMENSION IBUF(IOOOO) 

C COMPLEX AA ( 2 500 ) 

C CALL PLOTS ( I BUF ( 1 ), 1 0000.2 » 

C I PLOT * 0 

C PRD * 20.E-6 

DT = . 5E-9 
PI = 3.1415927 
F c 1.024E+6 


1203 

1204 

1207 

1500 

3007 

1502 

999 

1206 
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TF - 1.0E-6 
IPRINT * 10 
T * -DT 
WR J TE ( 6 * 1 1 ) 

11 format ( i hi j 
WRITE (6.12) 

12 FORMAT ( 10X.21HNUMERICAL INTEGRATION) 

WRITE <6*10) 

WR I TE ( 6 .2 1 ) 

21 FORMAT (6X.4HTIME.0X.5HlNPUT.8X*4HX( 1 ) *9x*4HX( 2) *9X*4 hX(3)*9X *4HX ( 4 
1 ) ) 

WRITE (6.22 ) 

22 FORMAT (6X,4HX< 5) .7X.4HX (6) .9X.4HX(7) *9X.4HX<8 ) ) 

WRI TE l 6.8 ) 

7 CONTINUE 

T * T + DT 

IPRINT * IPRINT ♦ 1 

U C 1 ) * 1. 

DO 4 I » 1 . N 
2D0T ( I ) * 0. 

DO b J = l.N 

5 ZDOT(I) * ZDOT(I) ♦ A(I.J)*2CJ) 

DO 4 J * 1»M 

4 ZDOT(I) = ZDOT(I) + 9 ( I • J ) *U ( J ) 

DO 6 I s l.N 

6 Z ( I ) a Z ( I ) + DT*ZDOT ( I ) 

I F t IPRINT. LE. 19) GO TO 7 
IPRINT * 0 

DO 2 I « 1. NORIG 
XII) * 0. 

DO 1 J a l.N 

1 xm ■ xii) + c u * j » *z i j ) 

DO 2 J a l.M 

2 XII) - XII) + Dll . J )*U( J > 

WRITE (6*8) 

WRITE ( 6.8) T«U(1).(X(I) * I *1 .NCRIG) 

C I PLOT = I PLOT + 1 

C. AA(IPLOT) = CMPLX t X ( 1 ) *0. ) 

8 FORMAT ( ) 

10 FORMAT ( 1H ) 

IFCT.LE.TF) GO TO 7 
C CALL TIMPLT (AA.IPLOT..PRD) 

999 CONTINUE 
RETURN 
END 

SUBROUTINE GREG 

COMMON A(13. 13)*B(13, 13). NOR IG.M*NfX(13)*NEVTeR.C(13. 13) *0(13*13). 
1 L AMREL.LAMIMG.se ALE *PV( 13 I .NLESS1 . ITEMP*Z (13) . I ROOT , I PRTFG 
REAL LAMREL .LAM IMG 
DOUBLE PRECISION PV 

C • t ******* t * « t # * * t *•»*»» t *«*»«»• t » 

REAL MATTMP ( 13 *13 ) 

DIMENSION CTEMPI 13.13 ) .DTEMP C 13 . 1 3 ) 

XLAM = 1. /LAMREL 

C ***** CALCULATION OF CTEMP ****** 

DO 2101 I * 1.NLESS1 
DO 2100 J = 1.NLESS1 

2100 CT£MP( 1 * J ) = -XLAM*PV( I )*A(N* J)*SCALE 

2101 C TEMP (1*1) = CTEMPCI.I) + 1. 

DO 2102 J = 1 « NLESS 1 

2102 CTEMP(N.J) = -XLAM* A ( N . J ) *SC ALE 
DO 2103 J a 1.NLESS1 

TEMP = CTEMP(N.J) 

CTEMP(N.J) = CTEMP( ITEMP. J) 

2103 CTEMP( ITEMP ,J) = TEMP 
6 FORMAT ( 1H ) 

2 FORMAT ( ) 

c ***** * *CALCUL AT I ON OF DTEMP 
DO 2110 I a l.NLESSl 
DO 2110 J = l.M 

2110 DTEMP ( I * J ) a -XLAM*PV( I ) *8 ( N* J ) 

DO 2111 J * l.M 

2111 DTEMP ( N ♦ J ) = -XLAM*B(N.J) 

DO 2113 J = l.M 

TEMP * DTEMP(N.J) 

DTEMP(N.J) * DTEMP ( I TEMP* J ) 

2113 DTEMP ( ITEMP, J) = TEMP 
COMMENT D * D ♦ C*DTEMP * * * * 

C C « C*C TEMP * * * * 

C SEQUENCE OF THESE TWO STEPS IS CRITICAL 

C * * * * UPDATE OF MATRIX D * * * * 

DO 3000 I * 1, NORIG 
DO 3000 J * 1 *M 
DO 3000 K « 1 .N 

3000 D(I.J) = D ( I , J ) ♦ C ( I .K )*DTCMP(K,J) 
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C * * * * UPPATF or MATRIX C ************ 

DO 3 001 I * I.NGRKj 
DO 3001 J - 1 ,N 
3UU1 MATTMPf I » J ) s C ( I * J ) 

DO 3002 I * l.NGRIG 
DO 3002 J = 1 « NLE S Gl 
ClhJ) = 0. 

DO 3002 < = 1 • N 

3002 ■ C ( 1 * J ) * C(I.J) * VATTMPI I ,<)*CTEMP(K.J) 

IF ( IPRTFG.EQ.O ) GO TO 8C08 
WRITE (6*2) 

WRITE (6* 100) I ROOT 

100 F0RVATU0X*9HMATRIX C ( * I 2 * 1H ) I 
DO 6502 I = 1 » NOR I G 

6602 WR I T E ( 6 * 2 ) (C(I.J).J * 1.NLE5S1) 

WRI TE( 6 *2 ) 

WR I TE ( 6 * 10 1 ) IROOT 

101 FORMAT ( 10X .9HMATRI X D(.I2*1H)) 

DO 6503 I = 1 .NORIG 

6503 wR I T E ( 6 * 2 ) (0(I*J)*J * 1*M) 

WR I T£ ( 6 * 2 ) 

8008 CONTINUE 
RETURN 
END 

SUBROUTINE DATAIN 

COMMON Afl3»13)*B (13*13) * NOR I G * v * N * X ( 13 ) .NEVTSR.Cl 13,13) *D (13*13)* 
1 LAMREL,LAMIMG,SCALE*PV( 13) * NL ESS1 * I TEMP.Z (13) . I ROOT * I PRTFG 
REAL LAMREL *LAMIMG 
DOUBLE PRECISION PV 

C*************** ****** ** ***•*•**«>* 

REAL LFACT1 . LFACT 2 . LT 1 . LV » LT2 , L 3 *L4 . 15 *L2 
IPRTFG = 0 
NEVTBR = 2 
N * 13 
M * 1 

C ***** INPUT PARAMETERS ********* ******** 

R 1 = 51. 

L T 1 = -5.358E-6 
LM = 6.06E-6 
LT2 = 5.194E-5 
RT2 * 9.11 
R2 = 5.1E3 
Cl = 4.8E-10 
L2 = 1.08E-5 
RL2 = 1.7 
C2 = 1 • 36 e £-9 
C 3 = 5.1E-10 
L3 = 5.1E-5 
RL3 = 8.01 
LA = 1.6E-5 
RL4 = 2.51 
C4 = 2.2E-9 
C5 = 3.44E-10 
L5 = 5.1E-5 
RL5 = 8.0] 

C6 = 5.7E-10 
R3 = 4.7E3 
C7 = 1.8E-9 
RL = 6.2E2 

LFACT2 = 1 ./(LT1*LT2*L v# (lT1+LT2) ) 

LFACT1 = LFACT2*(LT2*Lf-M /LM 
WRI TE (6*1500) 

1600 FORMAT ( 1 H 1 ) 

C * * * * * LOAD MATRIX A************ 

DO 1091 I = 1 « N 
DO 1091 J = 1 , N 
1091 A( I * J ) = 0. 

A(l. 10) - 1 ./Cl 
A( 1 .11) = 1 ./Cl 
A ( 2 * 1 0 ) = 1./C2 
A(3 .11 ) = 1 ./C3 
A(4 .11 ) = 1 ./C4 
A ( 4 * 1 3 ) = -1 . /C4 
A ( 5 . 1 1) = 1 ./C5 
A ( 5 * 1 2 > = -1./C5 
A ( 5 * 1 3 ) = -1./C5 
A ( 6.6 ) = -1 ./ { R3*C6 ) 

A ( 6 * 1 3 ) = 1./C6 
A ( 7. 7 ) = -1 ./ ( RL*C7 ) 

A( 7. 13 ) = 1 ,/Cl 
A ( 8 . 0 ) = -LM*R1 *LF ACT 1 

A ( 8 * 9 ) = -(RT2+R2)*(LT1«‘L* , )*LFACT1 + CR2*RT2)/LM 
A { 8 * 1 0 ) = R2* ( LT1*LM ) *LFACTl-k2/LV 
AI8. 11) = R2* ( L T 1*LV ) *LF AC T 1 -R2/LM 
A ( 9 * 8 ) = -LM*R1 *LFACT? 

A ( 9 ♦ 9 ) * -(RT2+R2 ) * ( LT 1*1.m ) *LFACT2 
A ( 9 , 1 0 j = R?* { L T dLM ) *L r AC T2 
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AI9.U) ■ R2* ( LT1+LM ) *LFACT2 

AI10*1) * -1./L2 

A ( 1 0 * 2 ) = -1./L2 

A ( 10*9 ) * R2/L2 

A ( 10*10) = -{ R2+RL2 I /L2 

A(1G,U) = -R2/L2 

Alllil J * -1./L3 

A ( 1 1 ♦ 3 ) = -1./L3 

At 11*4) = -1 • /L3 

A ( 1 1 * 5 ) = -1./L3 

A ( 1 1 * 9 ) * R2/L3 

A ( 1 1 * 1 0 ) * -R2/L3 

A ( 1 1 1 1 1 ) « -<R2*RL3)/L3 

A { 12 . 5 ) * l./LA 

A(12*12> * -RL A/LA 

At 1 3 * A ) * 1 • /L 5 

At 13*5 ) « 1./L5 

A ( 13 *6 ) *= -1./L5 

AC 13*7 | * -1./L5 

A C 1 3 • 1 3 > « -RL5/L5 

C * * * LOAD MATRIX B * * * * • «#••#« 

DO AOOO I * 1 • N 
DO AOOO J « 1,M 
AOOO BtlfJ) » 0. 

B ( 8 • 1 ) * LM*LFACT1 
8(9*1) * LM*LFACT2 

C * * * * **«# INITIAL CONDITIONS #*#***•** 

DO A001 I * 1 *N 
AG01 XU) - 0. 

C # * * * # PRINT INPUT DATA **#♦#** 

WR I T£ ( 6 * 100 ) 

100 FORMAT ( 10X.A8HTHE INPUT SYSTEM IS OF THE FORM XDOT • AX ♦ BU> 
WRITE (6. 6) 

WRITE (6*101) N*M 

101 FORMAT t 10X*28HT HE INPUT LINEAR SYSTEM HAS .I2.35M STATE VARIABLE C 
1CMPCNENTS AND HAS .I2.19H CONTROL COMPONENTS) 

wRI TE( 6.6) 

6 FORMAT (1H > 

2 FORMATt ) 

WR I TE ( 6 * 6 ) 

WR I T E ( 6 . 1 02 ) 

102 FORMAT ( 10X.3AHTHE INPUT MATRIX A IS AS FOLLOWS) 

WR I TE ( 6 *6 ) 

DO 200 I * l.N 

200 WR I T E ( 6 • 2 ) ( A ( I , J ) , J* 1 ,N ) 

WR I TE t 6 *6 ) 

WRITEJ6. 103 ) 

103 FORMAT ( 10X.3AHTHE INPUT MATRIX B IS AS FOLLOWS) 

WRI TEC 6 .6 ) 

DO 201 I * l.N 

201 WR I TE ( 6 * 2 ) <6(I,J)*-r • 1*M) 

WRITE (6*6) 

WR I T E ( 6 • 10A ) 

10A FORMAT ( 10X.A2HTHE INITIAL STATE VECTOR X IS AS FOLLOWS) 

WRITE (6*6) 

WRI TE ( 6*2 ) ( X ( I ) * I » l.N) 

WRITE (6.6) 

WR I TE ( 6 . 105 ) NEVT3R 

105 FORMAT ( 10X .A9HTHE SYSTEM IS TO BE SIMPLIFIED BY THE REMOVAL OF ,12 
1.12H EIGENVALUES) 

RETURN 

END 


149 



PRECEDING PAGE BLANK nos HLMffi 


APPENDIX I 

ROOT REMOVAL FOR A COMPLEX CONJUGATE PAIR 


151 



APPENDIX I 


ROOT REMOVAL FOR A COMPLEX CONJUGATE PAIR 


Although a complex conjugate pair of fast characteristic 
roots can be "removed" one at a time with ALSAP as described in 
the body of the report, the use of complex arithmetic in the 
simulator can be avoided by approximating the two modes simul- 
taneously. This appendix presents the equations describing the 
reduced state equation, the initial conditions for the new state 
variables, and the transformation relating the original and the 
new state variables. 

The original system will be of the form 


x = A x + B u , 

T 

where x = (x- , x 0 , ... x ) is the state vector and u = (u. , u n , 
1 ^ n l 2. 

... u^) is the control vector. Let the fast pair of complex 
conjugate roots be X and X, and let their eigenvectors be v and 
v, respectively (the over-bar designates the complex conjugate). 
As before, let the eigenvectors be normalized and partitioned 
such that 


v = 


2 1 


— 

and v = 

1 





1 i 


The reduced state equation will be 


.n 

x 


= A 


n n 
x 


+ B 


n 


u 


where the new state vector is given by x* = (x^, x^, ... > and 

the original state vector can be recovered using the relationship 

x = C x n + D u . 
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The initial conditions for the new state variables and the elements 
of the matrices A n , B n , C, and D are given by the following expressions. 
The superscript "re" refers to the real part of the complex number 
while "im" refers to the imaginary part; subscripts refer to matrix 
element indices. 

Initial Conditions 


n 

x. 

1 


= X. 


1 



for i = 1,2, . . .n-2 , 



re im re 

ViAvJ* p i J 


X 

n 


New State Distribution Matrix, A 


n 


n 

a . . = 
1J 


a u + 


K re , im 'N im re 1 

p n-l /p „-l > i - p i J 


n J 




^ > J 1,2, ...,n— 2. 


New Control Distribution Matrix, B n 



, im, im N . , T r im/ im N re 

b. . - l P. /p , )b . . + IP. /p , P t 

xj V i n- 1 s n-1 , j L V i r n-l S n-1 


- p rJ 


n J 


x 1,2 , . . . ,n-2 , j 1,2,..., m . 


State-to-State Relating Matrix, C 

The elements of the C-matrix are best described in three 
separate groups. 
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For i = 1,2 , . . . ,n-2 and j = l,2,...,m one has 

im re « re im 


C ij 1^2 { a n-l,j [ 


- JLUl J- U , 1C -L V 

X Pi - x p,. 


im 

? n-l 


re 


, T *n-l f ,re im , im re ^ Aim im 

+ 3 nj L — { X p i - X p i P i 


n-1 


+ X re P^ 6 } ] + | X | 2 6^ } , where 


'« - { 


0, if i + j 

1, if i = j. 


For i = n-1 and j = l,2,...,m one has 


, im re , re im 
X P„ X p„ 


= — 1 — J a P p n-r 1 p n-l "l. 
n-l,j 2 l n-l,j L im J 

A P n-1 


, im r im 
a . X P , 
nj L n-1 


(Ci ) <Ci ) ] } 


For i = n and j = 1,2, ••.,m one has 


1 J . re , . xm f re im \ 

: . = o 1 a . + X la ... - a .p )/ p f . 

nj i. \ z nj V n-1, j nj n-1 S n-1 J 


Control-to-State Relating Matrix, D 

The elements of the D-matrix are also described in three groups. 
For i = l,2,.*.,n-2, and j = l,2,...,m one has 
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For i 
d 

n 

For i 


, _ -1 f, r. re r re re im, im "\ 

d lt ' TTz t b ni V- V P 1 - Vl P i /p n-l ) 


Ul 2 1 

nj l 


' im 
s. p i 

+ \ ( 

+ b - 
n-1 

,[ 




- n-1 and j = l,2,...,m one has 

■ 1.1 - fl kj [»“ Vl ( [ p n-l ] 2 - Kl ] 2 ) ' (£l )] 


+ b . . [x re - x lm p re . / p 1 " ] } . 

n-l,j L n-1 n- 1 J J 


= n and j = 1,2,..., m one has 

d nj -^Kj [*“♦**“ to /*)] 
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